Brane world cosmology with Gauss-Bonnet and 

induced gravity terms 



by 

Richard A. Brown 



The thesis is submitted in partial fulfilment of the requirements for 

the award of the degree of 
Doctor of Philosophy 

OF THE 

University of Portsmouth 



November, 2006 



Copyright 

© Copyright 2006 by Richard A. Brown. All rights reserved. 

The copyright of this thesis rests with the Author. Copies (by any means) either in 
full, or of extracts, may not be made without the prior written consent from the Author. 



Declaration 



Whilst registered as a candidate for the above degree, I have not been registered for any 
other research award. The results and conclusions embodied in this thesis are the work 
of the named candidate and have not been submitted for any other academic award. 



To my family and friends. 



Edgar Allan Poe (1809-1849) ponders the subject of matter undergoing collapse: 

In sinking into Unity, it will sink at once into that 
Nothingness which, to all Finite Perception, Unity 
must be-into that Material Nihility from which alone 
we can conceive it to have been evoke d-to have been 
created by the Volition of God. 

from Eureka - A Prose Poem (1848), written 81 years before Edwin Hubble's observa- 
tions of the expansion of the universe. 



Abstract 



In this thesis we investigate certain cosmological brane world models of the Randall- 
Sundrum type. The models are motivated by string theory but we focus on the phe- 
nomenology of the cosmology. 

Two models of specific interest are the Dvali-Gabadadze-Porrati (DGP, induced-gravity) 
model, where the brane action is modified, and the Gauss-Bonnet model where the bulk 
action is modified. Both of these modifications maybe motivated by string theory. 

We provide a brief review of Randall- Sundrum models and then consider the Kaluza- 
Klein modes on Minkowski and de Sitter branes, in both the two and one brane cases. The 
spectrum obtained for the de Sitter branes is a new result. We then consider a Friedmann- 
Robertson- Walker brane in order to investigate the cosmological dynamics on the brane. 

We present a brief discussion of the DGP and Gauss-Bonnet brane worlds. We then 
investigate the Gauss-Bonnet-Induced-Gravity (GBIG) model where the Gauss-Bonnet 
(GB) bulk term is combined with the induced-gravity (IG) brane term of the DGP model. 
We present a thorough investigation of cosmological dynamics, in particular focusing on 
GBIG models that behave like self-accelerating DGP models at late times but at early 
times show the remarkable feature of a finite-temperature Big Bang. We also discuss the 
constraints from observations, including ages and Big Bang nucleosynthesis. 
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Chapter 1 
Introduction 



Quantum mechanics and General Relativity are two very successful and well validated 
theories within their own domains. The problem is that we have no way of unifying them 
into a single consistent theory. One of the most promising models of unification is String 
Theory (for an introduction to some of the concepts of String Theory see [HI). String 
Theories remove the infinities that are present in a classical unification by describing 
particles as extended 1 -dimensional strings rather than point particles. These ID strings 
live in a 10 dimensional space, or 11D for supergravity. The idea of extra dimensions 
was first put forward by Kaluza [2] and Klein |[3l in order to unify general relativity and 
electromagnetism. A perturbative approach is usually applied when working in String 
Theory due to the complexity of the equations. It has been shown that there are five 
different formulations of perturbative String Theory and that these are dual to each other, 
under a certain set of transformations. It has been suggested that these five theories are 
aspects of one underlying theory, M Theory. 

Two classes of strings are the closed and open strings. Gravity is described by closed 
strings and matter is described by open strings. In non-pertubative string theory there 
exist extended objects known as D (Dirichlet) branes. These are surfaces where the open 
strings must start and finish. This provides an alternative to the Kaluza-Klein approach, 
where matter penetrates the extra dimensions, leading to strong constraints from collider 
physics. If matter is confined to a 3-dimensional brane, the extra dimensions can be 
larger, since the constraints on gravity are weaker. 

In this thesis we consider brane world models. These models are inspired by String 
Theory, in particular by the Hofava and Witten model J4J. Brane world models are char- 
acterized by the feature that standard model matter is confined to a 1+3 dimensional 
brane while gravity propagates in the higher dimensional bulk. This means that gravity is 
fundamentally a higher dimensional interaction and we only see the effective 4D theory 
on the brane. We can not see or measure (with existing technology) Planck scale compact 
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dimensions as they require Planck energy scales in order to probe them. But if an extra 
dimension is large (relative to the Planck scale) then there may be signatures of this in 
colliders and table-top experiments, as well as in cosmological observations. 

In the model of Arkani-Hamed, Dimopoulos and Dvali (ADD) Q, the weakness 
of gravity on scales > 1mm (Newton's law has only been verified down to a scale of 
~ 0. 1mm ) is due to the existence of n ^ 2 extra compact dimensions. This model can 
be embedded in String Theory and has the advantage of putting forward a mechanism 
for explaining the hierarchy problem. The hierarchy problem is the discrepancy between 
the gravity scale (M P i anck ~ 10 19 GeV) and the electro-weak scale (M E w ~ 10 3 GeV). 
The 10 16 orders between these two scales is too vast to be considered natural. In the ADD 
model, length scales smaller then the scale of the extra dimensions (r < L) have a 4 + n 
dimensional gravitational potential: 

V~ r - (1+n) . (1.1) 

On scales larger than L, r ~ L in the n extra dimensions. Therefore the potential is given 
by: 

V~L- n r-\ (1.2) 

on large scales. The observed Planck scale is then the product of the fundamental Planck 
scale and the volume of the extra dimensions [[El : 

M 2 Planck ~ M\X n n L n . (1.3) 

Therefore in the ADD model the fundamental Planck scale could be much smaller than 
the Planck scale and the observed value is due to the size of the extra dimensions. 

The Randall- Sundrum model of two, Z 2 symmetric, 3-branes living in an anti-de 
Sitter (AdS) spacetime provides an alternative to the ADD model for n = 1 [9] (see 
Fig. 11.11) . The second Randall- Sundrum model has only one brane [ITOl . 

The Randall-Sundrum 2 Brane (RSI) Model 

The RSI model has two Minkowski branes living in a 5D AdS bulk, Fig. 11.11 The bulk 
metric is given, in a Gaussian normal coordinate system, by: 



ds 2 = a 2 {y)r] piV dx il dx v + dy 2 , 



(1.4) 



3 



TeV "hidden" brane 





Planck "visible 






-X 




5* \T 





:" brane 



Figure 1.1: The RS 1 configuration 

where rj^ is the Minkowski metric. We use Latin indices for 5D spacetime and Greek 
indices for 4D spacetime throughout this work. a(y) is the warp factor, given by: 

a(y) = e-\*\l t , (1.5) 

where t is the AdS curvature scale. The extra dimension is compact. In order for the bulk 
to be compact we invoke Z 2 symmetry so that: 



U++-y, L + y ^ L-y, (1.6) 

where L is the inter-brane separation. 

As we shall see in the next chapter the branes have equal and opposite tensions. At 
y = is the positive tension, "TeV" or "hidden" brane. At y = L is the negative tension, 
"Planck" or "visible" brane. The brane tensions A are given by: 

6 

\=o = ~\=l = -27 = A, (1.7) 

where k\ = 8ttG 5 = M^ 3 . We reside on the negative tension brane in order to solve the 
hierarchy problem. The "hidden" brane has fundamental energy scale M 5 . Then due to 
the warping of the bulk the effective fundamental energy scale on the "visible" brane is 
M Planck , where: 



M Planck = M\i [1 - e- 2L / £ ] 



(1.8) 
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y 



Figure 1.2: The RS2 configuration 



The low energy effective theory on the branes is of Brans-Dicke type [1 1 J. The Brans- 
Dicke parameter (u>) on each brane has the same sign as the brane tension [12J. Nucle- 
osynthesis constrains to > 40000 [fT3l[T4l . This is a problem due to the fact that in this 
model we reside on the negative tension brane and as such we should measure a negative 
value of u. A radion is required to stabilise the inter-brane separation in order for 4D 
General Relativity to be achieved at low energies 03). 

The Randall-Sundrum 1 Brane (RS2) Model 

The one brane (RS2) |[T0l model has a single positive tension brane at y = as seen in 
Fig. 11.21 The bulk is infinite in extent. As there is only one brane there is no mechanism 
for solving the hierarchy problem but the low energy theory on the brane is general rela- 
tivity lfl2l . Gravity is localised to the brane, in the low energy limit, via the warping of 
the bulk. At high energies the warping is insufficient to localise the gravity to the brane 
and observers on the brane perceive the 5D nature of gravity. The fundamental energy 
and the Planck scales are related by: 



i.e. we have effectively sent the second brane off to infinity (L = oo). 

The Randall-Sundrum model exhibits its 5D nature, to a brane observer, at high en- 
ergies, when the warp factor is insufficient to localise gravity to the brane. The deviation 
between the 4D and 5D nature can be seen in the gravitational potential. The modified 
gravitational potential is lfT2l : 



Ml 



Planck 




(1.9) 




(1.10) 
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In Ref. [ 10] this calculation was carried out with the result having a factor of 1 instead of 
2/3 in front of the r~ 2 term. The result in lfT2l was obtained by including the coupling of 
the graviton to matter on the brane. 

At high energies we can not get away from the 5D nature of the universe but at low 
energies (r ^> £) we are given a constraint on the AdS curvature scale. As we have seen 
no deviation from Newton's laws down to rs 0.1mm in the lab, Eq. (11.101) gives us the 
constraint that i < 0.1mm. This constraint means that the brane tension, Eq. (|1.7I) . and 
fundamental energy scale, Eq. (11.91) . are also constrained. These constraints are: 



A > (1 TeV) 4 , M 5 > 10 5 TeV. 



(1.11) 



In the RS models the 4-D cosmological constant on the brane is given by: 



A 4 = - (A 5 + K4A) . (1.12) 
Since the branes are Minkowski, we require the fine tuning: 



A 5 + «£\ = 0, (1.13) 

where n\ = 871G4 = Mpf anck , so that A 4 = 0. Breaking this fine tuning leads to different 
cosmologies on the branes 0: 



A > 



Ak 



K 1 



Anti-de Sitter, 



Minkowski, 



de Sitter. 



(1.14) 
(1.15) 
(1.16) 



The bulk cosmological constant is given by: 



6 

A 5 = -J 2 - 



(1.17) 



where i is the AdS curvature scale. 

The Randall- Sundrum models provides a simple way of investigating aspects of the 
phenomenology of the very complicated String Theories. In this thesis, we focus on 
brane world cosmological models based on generalisations of the RS scenario (see the 
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reviews dEl H61 [TTl [T8l [191 |20l [2T1). These include models that can explain the late-time 
acceleration of the universe without the need for dark energy. 



Chapter 2 

RS Brane Worlds 



2.1 Introduction 

In this chapter we will look in some detail at the Randall-Sundrum (RS) model. We 
will start with the simplest model, the Minkowski brane. We shall then consider a de 
Sitter brane before looking at a cosmological solution (the Friedmann-Robertson- Walker 
brane). In the first two cases we will look at the equations in the bulk and on the brane. 
We will also look at the spectrum of Kaluza- Klein modes generated due to perturbations 
on the brane. For the cosmological brane we will only consider the equations on the 
brane and how they differ from the standard general relativity result. It is these brane 
Friedmann equations that will be of interest in the following chapters. 

2.2 Minkowski Branes 

Minkowski branes are the most simple branes possible. They are flat, empty and static. 
We start by obtaining the background form of the bulk Einstein field equations. Using 
the bulk field equations we can work out the induced equations on the brane. Later we 
shall consider bulk gravitational waves. These give rise to the Kaluza- Klein modes on 
the brane. 

2.2.1 The background 

The line element for Minkowski branes in an anti de Sitter bulk was considered in the 
previous chapter, Eq. (11.41) . For simplicity we take y > and implicity assume Z 2 
symmetry. We can write the metric for this model as: 



g ab = a 2 (y)rj,J^ b +Sy a S!. 



(2.1) 
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This is a solution of the 5 -dimensional Einstein equation: 

G { 2 = -AbsS?- (2-2) 
Using Eq. (11.51) we can write out the following relationships for the warp factor: 



(2.3) 



where a prime represents differentiation with respect to y. This implies: 



and similarly: 



By Eq. (TO : 



2 

9ab,y = g'ab = —„ {9ab ~ W) , (2.4) 



fl* = *2, =>■ 9 ay = K. (2.6) 



These will be useful later. 

In order to construct the Einstein field equations we require both the Ricci curvature 
tensor and the Ricci scalar. We start by obtaining the RS Riemann tensor. The only 
non-zero derivatives of g a f, are those with respect to y. Using Eq. (12.41 ) we can write the 
Christoffel symbols as: 



Ite = - 7 (« + W - %9bc - W) , (2-7) 



so that: 



C = -|W(^-W)- (2-8) 
The Riemann tensor is defined as: 

r>a> pa pa i "ne pa ne pa q\ 

n bcd — 1 6d,c 1 6c,d 1 fed 1 ec 1 be 1 ed' 

So for the Randall-Sundrum brane model the Riemann tensor has the anti de Sitter form: 

R a bcd = -\ 2 (S a c g bd - s a d g bc ) ■ (2.10) 
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The RS Ricci tensor and scalar are given by: 

4 

Rab = —-jp9ab-, (2.11) 

and: 

R = -™- (2-12) 

The Einstein tensor is given by: 

(*ab = ~-^9ab ~ ~ 9ab = ~^9ab, (2.13) 



which satisfies the field equation: 

G ab = -A 5 g ab , (2.14) 
where the 5D bulk cosmological constant is given by: 

A6 = -|. (2.15) 
This also means that we can write the Ricci tensor as: 

2 

Rab = ^59ab- (2.16) 

Using the solution for the RS bulk cosmological constant in Eqs. (11.141) we find that 
in order for the brane(s) to have a Minkowski geometry we require the brane tension(s) 
to be fined tuned to: 

A± = ± ^|' (2 - 17) 

where the ± refer to the two different branes. The brane tensions can also be found from 
the extrinsic curvature of the brane. The extrinsic curvature depends on how the brane 
is embedded in the higher dimensional space. As we are working in Gaussian normal 
coordinates we can write the extrinsic curvature as 11221 : 

= \d y g^. (2.18) 
For Minkowski branes we find that the extrinsic curvature is given by: 

K^ = ^f- (2.19) 

The extrinsic curvature on the brane can also be obtained by considering the Israel- 
Darmois junction conditions 11231 : 
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K 



T, 



brane 
av 



-T 



brane 



9nv 



(2.20) 
(2.21) 



where: 



rpbrane rp 

1 (j.v — 1 uv 



(2.22) 



T^ u ane is the total energy-momentum tensor of the brane. This is the combination of the 



standard matter T^ v on the brane and the brane tension A which acts as a brane cosmo- 



logical constant. 

Using Z 2 symmetry, K~ v = ■ 
(where we drop the ± superscript): 



the extrinsic curvature can be written as O 



K. 



— -K* 



(2.23) 



where T = T£. Using Eq. (T2T231) for a Minkowski brane i.e T nu = and Eq. (I2TT8T) we 
find that the brane tension is given by: 



Ah 



±- 



6 



(2.24) 



Using this form of the brane tension with that in Eq. (12.171) gives us the following rela- 
tionship for Minkowski branes: 



(2.25) 



2.2.2 Field equations on the brane 



For a general cosmological brane, in a general bulk, the 5-dimensional field equations 
are: 



(5) _ k .2 T (5) 
ab — K '5 1 ab J 



(2.26) 



where the energy-momentum tensor is given by: 



rpa (5) rpa I _i_ rpa I 

J b -'ft \bulk ' -i b \brane- 



(2.27) 
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(5) 

The bulk part of T^ b is given by: 



A 5 

Tb\buik = nS b , (2.28) 

which is just the 5D bulk cosmological constant. 

The induced field equations on the brane are given by j8l : 

= -A 4 g^ + k\T^ + _ + 4y (2.29) 

The first two terms on the right hand side are the standard GR contributions of the matter 
on the brane and the cosmological constant. A 4 and k 4 are as defined in Chapter Q] The 
other terms are unique to brane models. The tensor is the high-energy correction 
term and is quadratic in the energy-momentum tensor: 

S, v = ^TT^ - -T^T? + 1^ [3T a pT a ? - T 2 } . (2.30) 

The fourth term in the induced field equations, S^, is the projection of the bulk Weyl 
tensor orthogonal to n a : 

£r = c2Ln e n i 9®%9 {B) l (2-31) 
The Weyl tensor describes the free gravitational field. The F^ v tensor is defined as: 



ab ^ 



9nu, (2.32) 



and includes contributions on the brane of the 5D energy-momentum tensor (if there is a 
bulk scalar field this is where it is felt). In the special case of Minkowski branes and an 
empty bulk, the Weyl tensor is zero, so that £^ u = and due to the empty bulk F^ v = 0. 
Also in Minkowski spacetime we have = and therefore = 0, and A 4 = 0. The 
induced field equations on the brane are simply given by: 

= 0. (2.33) 

2.2.3 Kaluza-Klein Modes 



In this section we shall look at the Kaluza-Klein (KK) modes on the Minkowski brane. 
Before we look at the solutions we should consider the polarisations of the 5D graviton 
as this generalises the equivalent 4D polarisations. 
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Polarisations of the 5D graviton 

The number of polarisations of the graviton is given by the number of degrees of freedom 
in the perturbation (/i a &) of the bulk metric tensor: 

9S "> 9% + K b . (2.34) 

The perturbation initially has 15 components since the metric is symmetric. We then 
apply the Randall-Sundrum gauge, which is the Gaussian normal and the transverse- 
traceless gauges combined. The Gaussian normal gauge (h ay = 0) removes 5 degrees 
of freedom. The transverse-traceless gauge (after the application of the Gaussian normal 
gauge) is: 

This removes another 5 degrees of freedom. This means that the 5D graviton has 5 
degrees of freedom and therefore 5 polarisation states. These polarisation states (as felt 
on the brane) are decomposed into [H: 

• 2 polarisations in 4D spin-2 graviton modes. 

• 2 polarisations in 4D spin-1 gravi- vector modes. 

• 1 polarisation in 4D spin-0 gravi-scalar modes. 

The normal 4D graviton is the massless spin-2 mode. The massive modes of the graviton 
are felt as massive modes in all of these fields on the brane, i.e massive graviton, gravi- 
vector and gravi-scalar. 

Obtaining the gravitational wave equation 

The perturbed line element is: 

ds 2 = (a 2 (y)^ + sh^v) dx^dx u + dy 2 , (2.36) 
where |e| <C 1. The modified metric can be written as: 

9ab = 9ab + ^9ab, 

= ga + eh^Sffi, (2.37) 
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where g ab is the background metric (all subsequent barred terms are background quanti- 
ties). Therefore for the inverse metric we have: 



g {5)ab = g ab -eh^6 a Jl 



(2.38) 



The perturbed Christoffel symbol is: 



£p(5)a 



6c 



■x {T d ( Vc + h dc , b - h bc4 ) - h ad (g bdiC + g dc>b - g bcA ) } . (2.39) 



The background Christoffel symbol is given by Eq. (12.71) , h ay = and g b = 0. This 
becomes: 



^ (5) I = ^ {f v (h v ,c + h uc , b - h bc>v ) - g ay h bc , y h™ (g bUjC + g^ b - g bCjU ) } . (2.40) 



Then by using the relationships in Eqs. (12.41) . (12.51) and defining h^ u = a 2 h^ u this can 
be written as: 



6T 



(5) a 
be 



K',c + h\, b - V - - lsfh,;< + sfhs 



c b ' b c 



: a 2 *: { l~h bc - hcv) + ? ( WV + W^a 



(2.41) 



In order to write out the perturbed Riemann tensor we need the explicit form of the 
derivative of Eq. (|2.41l) . We find this to be: 



ST 



(5) a 



bc,d 



h 7 4- h 1 
a b ,cd + a c,\ 



2r 



+ a 5y ( -h bC) d — h bCt yd 



2a 2 (2~ 

-k bc - h bC)V 



+ - ( W» a M + 



(2.42) 



We now have all the required quantities for the perturbed Riemann tensor which is 
given by: 



5R 



(5) a 



bed 



c r (5)a _ep(5)a =e ( 5 ) a =a (5) e -e (5 ) a -a (5) e 



bd,c 



bc,d ' bd ec ' ec bd be ed ed be' 



(2.43) 
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In order for us to obtain the perturbed Einstein equations we do not need to explicitly 
write out the full perturbed Riemann tensor. We require the perturbed Ricci tensor so we 
can just consider each term and perform the Ricci contraction on each in turn. 

We must also consider the Randall-Sundrum gauge requirements. We have already 
used h ay = but we also have h v „ v = and M = 0. Using these and the Ricci contrac- 
tion over the 1st and 3rd indices we find that Eq. (12.421) (the 2nd term in the perturbed 
Riemann tensor) has zero contribution to the Ricci tensor. Writing out Eq. (12.421) for the 
1st term in the perturbed Riemann tensor and carrying out the same procedure we find 
the perturbed Ricci terms to be: 



2 



: f ■ - 4a 2 r 4a 2 r 



~ V\ _ ~W hbd + T V» - a2h bd, yy \ • (2.44) 



The 3rd term in the perturbed Riemann tensor (T e bd 5T^ a ec ) also gives us zero when we 
perform the Ricci contraction. The 4th term in the perturbed Riemann tensor (rl c 5T^ e bd ) 
gives us the following components of the Ricci tensor: 

e f 8a 2 T 4a 2 7 1 
~ 2 I ~P bd ~ ~T M ' y J ' ( } 

The 5th term in the Riemann tensor (rg^r^ e bc ) gives us: 

s(2a 2 ~ 2a 2 ~ \ „ ^ 

~ 2 i I 2 M ~ ~T ' y J ' ( ^ 

The final term in the Riemann tensor (r^P ; ed ) also gives us the following terms in the 
Ricci tensor: 

e(2a 2 ~ 2a 2 , 1 

Then by combining Eqs. (12.441) . (12.451) . (12.461) . (12.471) and making use of the fact that we 
can change the indices as long as the summations stay correct we find the perturbed Ricci 
tensor to be: 

5R< 2 = I - -^r~ h bd + ^rhd, v - a 2 hbd, yy \ ■ (2.48) 



Using Eq. (12.1 II) we can write the following 



4 



SRf] = -^h d . (2.49) 
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By combing these two forms of the Ricci tensor we find that the following wave equation 
must be satisfied: 



M, 7 £ 

We can write the wave equation (12.501) as: 



An 2 - 



(2.50) 



ar*U + &,--d v 



= 0, 



(2.51) 



where □ is the 4D Minkowski d'Alembertian operator. This can be equivalently written 
with h bd instead of h^. The wave equation (12.511) is also: 



□ (5) V = o- 



(2.52) 



Using the wave equation (12.511) becomes: 



h^v = 0, 



(2.53) 



which agrees with the form given in j[2T|. The junction conditions can be enforced via a 
delta term which takes into account the presence of the brane at y = 0: 



hf„, = 0. 



(2.54) 



The two forms of the line element seen so far are written in Gaussian normal coordinates 
(the extra dimension is normal to a given hypersurface, the Minkowski spacetime brane 
in this case). Another coordinate system that can be used is the Poincare system. This 
gives us a line element of the form: 



ds 2 = — [dz 2 + rj^dx^dx^ 



(2.55) 



where z is defined as: 



(2.56) 



The perturbed line element in these coordinates is given by: 



ds 1 



dz + ( r]^ + sh^v ) dx^dx 



(2.57) 
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Solving the wave equation 

The solutions to the wave equation (12.531) can be written as the superposition of Fourier 
modes of the form: 

h^(x\y) =U n (y)e^ xX e MV , (2.58) 

where is a (constant) polarisation tensor. We are interested in U m (y), the bulk gravi- 
tational wave part of the decomposition. By using Eq. (12.581) in Eq. (12.511) we find the y 
dependent wave equation to be: 

e- 2y/e U m = 0, (2.59) 

where m can be interpreted as the 4D mass of the 5D graviton. The mass is due to the 
fact that an observer on the brane would measure the timelike 4D projection on the brane 
of the gravitons null 5-momentum: 

-m 2 = k x k x . (2.60) 

The zero mode satisfies the boundary conditions given by U = C where C is a constant. 
For m > 0, we now define the function: 

F m (y) = e- 2y / £ u m (y), (2-61) 

and change variable to: 

X = e y/e m£, (2.62) 

to rewrite Eq. (12.591) as: 

X 2 ^ + X d -^+(X 2 -A)F m = 0. (2.63) 
This is a Bessel equation, with solution: 

F m = Z 2 (e^me) , (2.64) 

where: 



Z P (S)=AJ P (S) + BY P (S), 



(2.65) 




Figure 2. 1 : A plot of the potential in equation (12.701) with t = 1 



with A and B constant and p is the order of the Bessel function. By Eq. (12.581) the full 
solution for m > is: 

K„(x x , y) = e 2 y' l Z 2 (e^W) e^e^. (2.66) 
Schrodinger form of the wave equation 

The wave equation (12.591) can be rewritten in a form similar to the Schrodinger wave 
equation. We define two new variables as: 

V> m = a 3 / 2 U m , (2.67) 

Z = £(e y/i - 1). (2.68) 



This gives: 



°~ l +V(Z)ip m = m 2 iP m . (2.69) 



dZ 2 

This is a Schrodinger equation with a potential V(Z) given by: 



1 ^ 

= WTW (2 ' 70) 



This is plotted in Fig. I2.ll 
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Figure 2.2: A plot of equation (12775b with L/£ = 1/4, 1/2 and Ljl = 1 for the solid, 
dashed and pale solid lines respectively. 



Specific Solutions 

In order to obtain specific solutions we must apply boundary conditions to the general 
solution. The boundary conditions for the perturbations are: 



dy 



0. 



brane 



(2.71) 



This means RSI has two boundary conditions whilst RS2 only has one. For RSI: 



A ^y ( e2 ^ J 2(e^W)) \ {y=Q) + B^- (e 2 y/%(e^ e m£)) \ {y=0) = 0, (2.72) 



A Yy (e 2v/e Me y/e mi)) \ {y=L) + s| (e^Y 2 (e^mi)) \ {y=L) = 0. 
These can then be combined in the matrix equation: 



(2.73) 



J 2 (0) K 2 '(0) " 




' A ' 


J' 2 (L) Y>(L) 




B 



0. 



where J 2 (0) = ^ ^(e^mf) |( y=0 ) and so on. 

In order for us to obtain a unique and non-trivial solution we require that: 



(2.74) 



4(o)y 2 '(L) - J 2 (L)y 2 '(o) = o, 



(2.75) 
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i.e.: 

Ji (mi) Y l (e L/e mi) - J x (e L/e mi) Y x [ml) = 0. (2.76) 

In Fig. 12.21 we see that for a particular value of L/i we have a series of values of mi 
that satisfy the above equation. We see that as we increase the brane separation the mode 
spacings decrease. We find that we have a zero mode solution with an infinite tower 
of discrete massive Kaluza-Klein modes. In the RS2 scenario we only have a single 
brane and therefore only one boundary condition. This means we can not solve for m 
independently of A and B and therefore obtain a specific solution. Thus RS2 has an 
infinite continuous tower of Kaluza-Klein modes. This is confirmed by taking the RSI 
solution and sending the second brane off to infinity. 

2.3 de Sitter Brane Worlds 

The de Sitter brane is the simplest expanding solution. A de Sitter universe expands at a 
uniform rate given by the Hubble constant H. 

2.3.1 The background 

In Poincare coordinates the spacetime is described by IfTTTl : 

ds 2 = A(z) 2 (dz 2 + >y ia ,dx' i da?) , (2.77) 

where: 

^ u dx^dx u = -dt 2 + e 2Ht dx 2 . (2.78) 

and the warp factor is: 

A(z) 

On the visible (positive tension) brane, 
solution (A(z + ) = 1), so that: 



sinh (Hz) ^ ^ 

z = z + , we need to have the standard de Sitter 



z + = H-hinh- 1 (Hi) 



(2.80) 
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Let us consider the two brane case where the second brane is at z = Z- > z + . The 
extrinsic curvature in this coordinate system is given by: 



leading to: 



= - l -A- x d z g^ (2.81) 



= . cosh^) (2 _ g2) 



The extrinsic curvature as found from the junction conditions, Eq. (|2.23l) . is independent 
of the brane geometry. For de Sitter branes = —p±g^ v with p ± constant, unlike the 
Minkowski case, where = 0. The extrinsic curvature for a de Sitter brane can then 
be written as: 

= -ycr^. (2.83) 

where a — A + p is the total energy density on the brane. By Eq. (|2.82l) we obtain the 
following relationship: 

a± = ±^- n cosh(Hz±). (2.84) 

k 2 5 £ 

There is some degeneracy in the brane tension and the energy density on the brane. We 
can rewrite the above equation in terms of the brane tensions as: 

A± = ±-?-cosh(Hz ± ) - p±. (2.85) 
Now if we evaluate the metric on each of the two branes: 

• Positive tension brane: 

ds\+ = -dt 2 + e 2Ht dx 2 . (2.86) 

This is just the standard de Sitter line element as we have set the observer on this 
brane. 

• Negative tension brane: 



ds%_ = A 2 {z^)(-dt 2 + e 2Ht dx 2 ) 1 

= -df+e^-t-dx 2 , (2.87) 
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where we have defined the proper time on the negative tension brane as: 

t_ = A(z_)t + C, (2.88) 

with C constant. The Hubble rate on the negative brane is also modified by the warp 
factor: 

This means that the warp factor on the negative tension brane is less than that on the 
positive tension brane 11241 . 

The form of the Hubble rates given above means that the two are dependent only on 
the brane position (z). We can view the expansion of the brane as a consequence of the 
brane moving through the warped bulk. We can write the Hubble constant on each brane 
as: 

sinhi/i z+ 

H± = (2.90) 

where H + = H (the Hubble rate on the positive tension brane). Using this result in the 
brane tension, Eq. (12.841) . we find that the tensions are given, as functions of the Hubble 
rate, by: 

\± = ±-\-Jl + {ZH±y-p±. (2-91) 
High energy branes are characterised by Hi ^> 1, so that: 

A±«±-^# ± -p±. (2.92) 



On low energy {Hi <C 1) branes: 



\ ± ^±—-p ± . (2.93) 



This reduces to the Minkowski relation when p± = 0. 
In the de Sitter case: 

S,u = -^P 2 9»v (2.94) 

The bulk is still empty in this case so the other two terms (£^ u and F^ u ) are still zero. The 
induced field equations on the brane are then given by: 

K 2 

= -A 4 g^u - ^\p9y,v - 7^P 2 SV' (2-95) 
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where 



2 _ ^ 4 



6 

Using Eqs. (11.121) and (12.151) we can write the field equations on the brane as: 



G 



\2 4 

6 



6 



1 flV 



(2.96) 



(2.97) 



2.3.2 Kaluza-Klein modes 

In Gaussian normal coordinates, the metric for a de Sitter brane is: 

ds 2 = dy 2 + w 2 (y)'y f j il/ dx ll dx u , 
where warp factor is given by: 

w(y) = H£smh(y/£). 
If yx is the position of the visible brane then: 

1 



H 



hinh( yi /£)- 



The perturbed metric is: 



ds 2 = dy 2 + w 2 (y) (7^ + f^) dx^dx\ 



where: 



(2.98) 



(2.99) 



(2.100) 



(2.101) 



Tf^ = = Vf^ 
The perturbation can be split into Fourier modes: 



(2.102) 



fw (t, x, y) -> fm(y)ip m (t,k) (x) . 



(2.103) 



Using these in the perturbed field equations we obtain the following two wave equa- 
tions [|2l: 



and 



ip m + 3H?p m + 



k 

7m 



+ m 



= 0, 



(2.104) 



(2.105) 
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where m is the separation constant and again can be considered as the effective mass of 
the 4D graviton. Eq. (12.1041) becomes: 



cosh(y/£) 



f' m (y) + 



( m V 
\Hhinh{y/£) ) 



£swh(y/£) 

The zero-mode satisfying the boundary condition is given by: 

fo = constant. 
For the massive modes, m > 0, we rewrite Eq. (12.1061) as: 



f m (y) = o. 



(2.106) 



(2.107) 



-d y {sinh\y/£)d y ) + 



(- 



m \ ' 



sinh 4 (y/£) yK ^ \H£swh{y/£)J 

The wave equation for (12.1051) can be written as: 

[□ - m 2 ] i) m = 0. 
If we make the following substitutions: 



f m (y) = o. 



(2.108) 



(2.109) 



U = smh 3 / 2 (y/£)4>(y), Z = cosh(y/£), 



we can write Eq. (|2.108l) as: 

d 2 U 



dU 

[ - z) !zt- 2Z lz + 



1-Z< 



U = 0, 



(2.110) 



(2.111) 



where v = | and ji = i^J ((m/H) 2 — 9/4) = ip. The solutions to this equation are 
associated Legendre functions. The general solution is a linear combination of the two 
associated Legendre functions and Q£. We find the profile for the KK modes are 
given by: 



Mv) = sinh- 3 / 2 (y/£) \ A p P^ im/H)2 - 9/4 



(cosh (y/£)) 



B pQ 3/2 (C0Sh (?//£)) 



(2.112) 



The two branes, at yi and y 2 , are related by: 



y 2 - yi = L. 



(2.113) 
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Figure 2.3: Determining the KK mass eigenvalues for two de Sitter branes with Ljl = 
1/4, 1/2 and Lj£ = 1 for the solid, dashed and light solid lines respectively. The y-axis 
is the left4iand side of Eq. (12.1161) . H£ = 0.1 in this plot. 



We assume yi > 0, y 2 > y\- The boundary conditions are the same as in Eq. (I2.71I ). To 
obtain forms for the derivatives we use the following relationship [|26l : 

d ^f L = {Z 2 ~ 1)^ {nZPftZ) - 0*+ u)P^(Z)} , (2.114) 
which is valid for Re(Z) > 1. Q$(Z) obeys the same relations. Then: 



df m (y) / 3 / 2 + *\/(™/H) 2 -9/4 



dy 



£sinh 5/2 (y/£) 



A m Plf m,H)2 - m (cosh (y/i)) 



+ B m Q% {m/H)2 ~ 9/ \cosh(y/i)) 



(2.115) 



We then obtain the specific solution to the boundary conditions: 



Pl P /2 (cosh ( yi /£)) Q\ p /2 (cosh (y 2 /£)) - Pfc (cosh (y 2 /£)) Qf /2 (cosh ( yi /£)) = 0, 

(2.116) 
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mi 

15 



u 1 2^3 4 5 

1 

Figure 2.4: A plot of the first massive mode as a function of L/£ (solid line), and the mass- 
gap (dashed line). We see that as we increase L/t the first massive mode approaches 
3H/2. Hi = 1 in this plot. 



Vl arcsinhT^V ^ = ^ + ^. (2.117) 



where: 



t \Hl r 

The KK mass eigenvalues defined by this equation are illustrated in Fig. 12.31 As far as 
we are aware, these eigenvalues have not previously been found. 

By rewriting Eq. l2.108l in Schrodinger form, it is possible to show that for normalis- 
able modes E51 : 

m 2 >^H 2 . (2.118) 



i^/ (m/H) 2 — 9/4 

In the general solution, Eq. (|2.112l) , the associated Legendre function P 3 J 2 
diverges for large y if iyj (m/H) 2 — 9/4 is real i.e. for m 2 < \H 2 . There is a mass gap 
between the zero mode and the normalisable KK tower in both the one and two brane 
set-up. 

As with the Minkowski solution if we increase the separation of the branes the mode 
separation decreases (see Fig. 12.31) . In Fig. I2.4l we can see the spacing between the zero 
mode and the first massive mode as a function of L/L If we consider the RS2 model and 
send the second brane off to infinity we find that there is again a continuum of massive 
modes starting at m 2 > \H 2 (m 2 = f if 2 is not a solution). 
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Figure 2.5: The KK mass eigenvalues for two de Sitter branes with different values of 
H£, i.e. different energy regimes. H£ = 0.1, 1 and Hi = 10 for the solid, dashed and 
light solid lines respectively. The y-axis is the left-hand side of Eq. (12.1161) . Ljl — 1 in 
this plot. 

In Fig. I2.5l we see the effect of going from the low to the high energy regime. At high 
energy the first KK modes are much more massive than in the low energy regime. The 
spacings are also much larger. Fig. 12.61 shows the spacing between the first massive mode 
and the mass gap as a function of Hi. As we increase the energy (Hi increasing) the 
mass gap increases linearly. The mass of the first KK mode also increases linearly but at 
a faster rate than the mass gap. Therefore at higher energies we must put in even more 
energy, relative to the low energy regime, in order to excite the KK modes. 

2.4 Friedmann Equations 

Here we shall consider Friedmann-Robertson- Walker (FRW) branes. This is a general 
cosmological solution where the Hubble rate decreases as the energy density goes down. 
The FRW solution describes the uniform expansion of a homogenous and isotropic per- 
fect fluid. 

We shall derive the Friedmann equations, that govern the evolution of the model, for 
the standard GR case before we consider the brane results. 
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Figure 2.6: A plot of the first massive mode as a function of Hi (solid line), and the 
mass-gap (dashed line). We see that as we increase Hi the first massive mode diverges 
from 3H/2. L/£ = lm this plot. 



2.4.1 The GR result 

The FRW metric is: 



ds 2 = -dt 2 + a 2 (t) 



dr 2 



+ r 2 d6 2 + r 2 sin 2 9 \ , 



(2.119) 



\-Kr 2 

where t = constant are maximally symmetric 3-spaces. The scale factor is a(t) and K 
is the spatial curvature. We then use the Einstein Field equations: 



For A = 0, k 2 = 8nG: 



(2.120) 



We use the energy-momentum tensor for a perfect fluid: 



(2.121) 



where is the 4-velocity: 



= (p + P) U»U V + Pg^, 



ds 



(2.122) 



(2.123) 
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Therefore: 

T = -p + 3P. (2.124) 
The spacetime is expanding but the fluid is at rest in a comoving frame i.e: 

[^ = (1,0,0,0). (2.125) 
The other ingredient that we need is the conservation of energy equation: 

V M T^ = 0. (2.126) 

The time component gives: 

p=-3pH(l+w), (2.127) 
where w = P/p is the equation of state (\w\ < 1). If w is constant then: 

poca- 3{1+w) . (2.128) 

Throughout this work we assume that w is constant. Physically this is not true as the 
effective equation of state will decrease from iy = l/3tow = 0as the universe evolves 
from radiation domination to matter domination. But is reasonable to assume w is con- 
stant as we either consider early time, radiation domination, or late time, matter domina- 
tion. 

By Eq. (12.1211) we find the pure temporal component of the Ricci tensor to be: 

#oo = 4vrG(p + 3P). (2.129) 
From the metric we can show that: 



i?oo = -3-. (2.130) 

a 



So we have: 



- = -^G(p + 3P), (2.131) 
a 3 

which is one of the Friedmann equations. We now consider the spatial directions. In fact 
we only need to look at one as the expansion is isotropic. If we look at the 99 component, 
from the metric we have: 



R ee = r 2 (aa + 2d 2 + 2K). 



(2.132) 
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R ij = 4nG(p-P)g lJ , (2.133) 

so the 99 component is: 

R ee = 4nG(p- P)a 2 r 2 . (2.134) 
By Eqs. (12.1321) and (12.1341) we have: 

- + 2^ + 2^ = 47rG(p-P). (2.135) 
a cr ar 

Using this with Eq. (|2.131l) we get the main Friedmann equation: 

a\ 2 8nG K 



p-- . (2.136) 



a J 3 a 2 
This is usually written in terms of the Hubble parameter which is defined as: 

H = -. (2.137) 

a 

If we include a cosmological constant in the field equations the Friedmann equations 
become: 

- = -4rG(p + 3P) + ^, (2.138) 
a 3 3 



and: 



rr2 8ttG K A 

# 2 = — + (2.139) 

3 a 2 3 



These are the standard GR Friedmann equations. When we consider the braneworld we 
will require the Israel junction conditions. This is because we now live on a hypersurface 
so the Friedmann equations only apply to a restricted part of the spacetime. 

2.4.2 RS result 

We consider a single Friedmann brane model (i.e RS2), so that gravity is described by 
General Relativity in the low energy limit (and not a Brans-Dicke theory). The line 
element is E71 : 

dsf 5) = -N(t, y)dt 2 + A 2 (t, y)-f ij dx l dx j + B 2 (t, y)dy 2 , (2.140) 
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where 7^ is the maximally symmetric space given by: 

dr 2 



^dx'dx 3 = + r 2 d6 2 + r l sin^ 9d^ '. (2.141) 

1 — Kr l 



The 5D field equations are: 



G£ (2.142) 



where the energy-momentum tensor is given by: 



T b — T b \bulk + T b a \ brane- 

(2.143) 



(5) 

The bulk part of T afe is given by: 



T b a \ bulk = -^51 (2.144) 
which is just the 5D bulk cosmological constant. The brane term is given by: 

T b a \ brane = ^diag(-p - A, P - A, P - A, P - A, 0), (2.145) 


where A is the brane tension. The standard conservation of energy equation applies on 
the brane (only gravity can leak off the brane into the bulk). If we follow [|27ll and look 
at the non-zero components of the Einstein tensor we get for the tt component (' = d/dy 
and' = d/dt): 



, 5) I A A B\ n 2 (A" A 1 (A' B'\\ N 2 . 



For the three brane dimensions we get: 



y(5 ) A 2 (A' (A' _N'\ B' fN' _A!\ _A" N" 



13 B 2 \ A \ A N J B \ N A J A N 

A 2 \ A ( N A\ A B ( N A\ B\ 
TP*\A \ 2 N - a) - 2 A + B [n - 2 A) - b) - -»' K - 

(2.147) 

For the bulk dimension we get two non-zero components: 

Gt y - 3 \a1v + ab~a^ (2 - 148) 
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and: 



Since = (and therefore = 0), this means that there is no flow of matter along 
the extra dimension. Therefore we get the following results from Eq. (12.1481) : 

B A! N' A 
In Ref. E71 they have shown that if we define the function: 

we can rewrite the Einstein equations in a simpler form. Using Eqs. (12.1461) . (12.1501) and 
(12.1511) we can show that: 

2A'A 3 

F' — — - — K 2 Tt\buik- (2.152) 
Using Eqs. (12.1491) . (12.1501) and (12.1511) we can show that: 



2AA 3 
~3 

Thus: 



P = ^^K 2 ry\ bulk . (2.153) 



■24V\buik f A , A3dy = _^As + c (2 154) 



3 J "6 
So if we now use the form of F from Eq. (|2.151l) we get 



2 -2 



AX /A ' Y K A5 + ^ (2.155) 



AN I \AB J A 2 6 A A 

We have the freedom in the coordinate system to redefine the extra dimension coordinate, 
so we set B = 1. We can also define t to be the proper time on the brane so we can set 
N(t, 0) = 1, which implies that: 
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We must now consider the A! /A term, and we use the Israel junction conditions: 



4 (- i 



K^ = —^\T llv --Tg llv ), (2.157) 
where: 

Kfu, = —=dy9ia>- (2.158) 
Using this we can show that we can write the extrinsic curvature of the brane as: 

Kit = diag (2.159) 



N A A A 



By Eq. (12T571) : 



A' 
1 



brane 



^(p + A). (2.160) 



By Eqs. (11.121) and (|2.17l) we can write the Friedmann equation as: 



where C is the constant of integration from Eq. (12.1541 ) (it can also be obtained from 
the electric (Coulomb) part of the bulk Weyl tensor [|28l ). This term is called the dark 
radiation term from the bulk due to its a~ 4 behaviour. It can be shown that C is the mass 
of the black hole in the Schwarzschild anti de Sitter bulk [|29ll30l . The main differences 
between this Friedmann equation and the GR result is that we have a p 2 term driving the 
expansion. When p ^> A, i.e. at early times (high energies), we see that H ~ p, unlike 
the GR case, H ~ v /p. Therefore the scale factor evolves like (taking C = 0): 

a oci 1/3(1+w) , (2.162) 

rather than: 



a oc 



t 2 '^ 1+w \ (2.163) 



as it does in the standard GR case. In the low energy limit p <^ \ (with C = again) the 
Friedmann equation gives us the same result as in the GR case. 

The Raychaudhuri equation for the FRW brane follows from the Friedmann equation 
(12.1611) and the conservation equation (12.1271) as: 

*=-f<" + »>( i+ &) + £-4- <2 - i64) 
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The RS form of the Friedmann and Raychaudhuri equations are going to modify the 
predictions of primordial nucleosynthesis due to the presence of the p 2 and dark radiation 
terms. In order for the RS model to mimic the GR result at nucleosynthesis we require: 

C 

A >> pnuc » -j. (2.165) 
We shall consider the nucleosynthesis constraint further in Chapter |3j 

2.5 Conclusions 

In this chapter we have looked at the field equations in the bulk and how they are projected 
onto the brane for Minkowski and de Sitter brane worlds. These equations could be much 
more complicated if we had considered a scalar field to be present in the bulk. In general 
a scalar field should be present in order to stabilise the brane positions when there are 
two branes. 

We have seen how the brane tensions are related on the two branes in both the 
Minkowski and de Sitter cases. The brane tensions are required to be fine tuned in order 
to obtain Minkowski or de Sitter branes. We have also seen how the Hubble rate on the 
two de Sitter branes are related. 

We have looked at the Kaluza- Klein modes on Minkowski and de Sitter branes in both 
the RSI and RS2 models. The Kaluza-Klein modes are decomposed into five different 
polarisations in three different fields (gravi-scalar, gravi- vector and the graviton). We 
have found the wave equation for bulk gravitational waves. We worked out the spectrum 
of the KK modes in both the Minkowski and de Sitter brane cases, the de Sitter case 
being a new result. In the RSI case there is a discrete spectrum, described by associated 
Legendre functions, beginning above a mass gap. The spectrum is altered due to the 
curvature scale and the inter-brane separation. In the RS2 model, with an infinite brane 
separation, we obtain a continuum of massive KK modes, starting above the mass gap. 
The mass-gap is a function of the Hubble rate on the brane. 

We then looked at a more general and physically relevant solution, that of the Friedmann- 
Robertson-Walker brane. We saw how the Friedmann equations were obtained in the 
standard GR setup and how we extend this to the brane world case. We saw that the 
brane setup modifies the high energy Hubble rate but at low energies we obtain the same 
result as in the GR case. In the next chapter we shall consider more complicated brane 
models which can give us more drastic modifications to GR. 



Chapter 3 

The DGP, GB and GBIG Brane World 
Cosmologies 

3.1 Introduction 

So far we have been considering the Randall-Sundrum brane models in which gravity 
appears 4D on the brane via the warping of the bulk dimension. This warping means that 
to a brane observer gravity appears 5D at high energy and 4D in the low energy regime. 
It is also possible to have 5D effects at low energy |j3Tll32l . The DGP (Dvali, Gabadadze 
and Porrati) model achieves this via a Einstein-Hilbert brane action; this brane Ricci 
scalar can be interpreted as arising from a quantum effect due to the interaction between 
the bulk gravitons and the matter on the brane ll32l (see also ll29l l33l ). The induced 
gravity term on the brane dominates at higher energies, below a certain length scale r, so 
that gravity becomes 4D at high energy. As the bulk is no longer required to be warped, 
the DGP model lives in a Minkowski bulk. When a Friedmann brane in used within 
the DGP model it is possible to show that one of the solutions gives rise to late time 
acceleration without the presence of a dark energy field li34l . 

We can also construct brane models using the Gauss-Bonnet (GB) higher order curva- 
ture terms in the bulk action (see, e.g., [[35l[lS[37l[M[32l4QllIQ^ 
These models modify gravity at high energy, unlike the DGP model. The Gauss-Bonnet 
term may be interpreted as a low-energy String Theory correction to the Einstein-Hilbert 
action. These models are akin to the Randall-Sundrum ones, in the sense that it is the 
nature of the bulk that gives rise to the 4D and 5D regimes on the brane, unlike the DGP 
model. 

We would like to have a model that has 5D effects at both high energy and low energy 
but 4D gravity in between. This is because we would like a model that can modify the 
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early universe and explain the acceleration of the universe that we are currently experi- 
encing. We require 4D gravity between these two limits in order for nucleosynthesis and 
other constraints to be obeyed. 

In this chapter we shall start by taking a brief look at the DGP and GB models and 
see how they modify the cosmology on the brane. We shall then consider the effect of 
combining these two. 



In 2000 Dvali, Gabadadze and Porrati (DGP) 11321 put forward a brane model that has 
low energy (infrared or IR) modifications through an induced gravity term. The induced 
gravity term (4D Ricci scalar term in the brane action) can be motivated by possible 
quantum effects of the interaction between matter on the brane and the bulk gravity. 
(String theories with a ghost free GB term in the bulk give rise to induced gravity terms 
on the brane [48].) 

At early times, the 4D term dominates and General Relativity is recovered (in the 
background - note that the perturbations are not General Relativistic [|49]l50ll5TI ). Grav- 
ity on the brane is 4D until the "cross-over" scale r is reached and we enter a 5D regime. 
At scales greater than r gravity "leaks" off the brane (the 5D Ricci term in Eq. (|3.1I) be- 
gins to dominate over the 4D Ricci term), and appears to observers on the brane to be 
5D. As we have 4D gravity on the brane from the induced gravity we do not require the 
warping of the bulk space in order to recover 4D gravity. So the DGP model lives in a 
infinite Minkowski bulk. The DGP model has two branches, one of which is very inter- 
esting for cosmology as it ends with a phase of "self-acceleration" [[341 . This is useful as 
we observe the universe to be in a period of acceleration which is usually explained via 
a dark energy field. The DGP model explains this with modified gravity. The DGP mod- 
els are in some sense "unbalanced", since they do not include ultra-violet modifications 
to cosmological dynamics. We would like a brane world model that modifies gravity at 
early times as this is where string and quantum effects must eventually dominate. 

With induced gravity on the brane we have a gravitational action of the form |[34l : 



3.2 DGP Branes 



s, 



grav 




(3.1) 
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where r is the cross-over scale and is defined as: 

4 Ml 

r = 4 = ^ (3 - 2) 

The sign of the induced gravity on the brane must be the same as the 5D part. We have 
an effective 4D gravitational constant which can be written as k\ = n\/2r. 

The bulk field equations are of the same form as in the RS model. The 4D energy- 
momentum tensor is modified, in comparison to the RS model, as it includes the con- 
tribution from the induced gravity term in the action. By using the same assumptions 
and conditions as in the RS case (a perfect fluid undergoing isotropic and homogenous 
expansion) we obtain the induced gravity Friedmann equation to be Il34l : 



where e = ±1 and C is the bulk black hole mass (the DGP limit has a Minkowski bulk 
A 5 = with C = 0). If we consider a Minkowski bulk (A 5 = 0, C = 0) with A = on a 
flat brane (K = 0) then we have: 

E 2 = ±- + | P . (3.4) 

r S 

The two different solutions corresponding to the two different values of e correspond 
to two different embeddings of the brane within the bulk. Both branches have a 4D limit 
at high energies: 

K 2 

DGP(±): H^r- 1 H 2 = ^-p, (3.5) 
while at low energies both have modified 5D limits: 

DGP(+): a -> H -> - , (3.6) 

r 

DGP(-): p^O H = 0. (3.7) 

DGP(-) has a non-standard (and non-accelerating) late universe. The self-accelerating 
DGP(+) branch is of most interest for cosmology. 
In dimensionless variables: 

h = Hr, ii = ^$p, r = -, (3.8) 
o r 

we can write the Friedmann equation (in the case C = A 5 = A = 0) as: 



h 2 = ±h + fj,. 



(3.9) 
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Figure 3.1: Solutions of the DGP Friedmann equation (h vs //) with C = A 5 = A = 0. 
(Brane proper time r flows from right to left, with r = oo at ji = 0.) 



The solutions for h can be written in terms of /x as: 



h=^±^± 1 -. (3.10) 
2 2 v 

These solutions are shown in Fig. 13. II 

In the + branch we have late time self-acceleration (H = 1/r) i.e we end in a vacuum 
de Sitter state. In the — branch we have H = so we end up with a Minkowski state. The 
DGP(+) model gives us late time acceleration without dark energy with early universe 
dynamics that agrees with GR so we do not need to modify baryogenesis and other early 
universe events. The DGP late time acceleration has been tested against the supernova 
luminosity and baryon acoustic oscillations Il52l[53ll54l . But the DGP is "unbalanced" as 
we would expect high energy (ultraviolet or UV) modifications as well. If we generalise 
the DGP model by introducing a negative bulk cosmological constant and a brane tension, 
so that the bulk is anti de Sitter, we do not get the desired results. It has been shown 11551 
491 [56ll that if we have an induced gravity brane in a warped bulk we either have 4D at 
early and late times with a period of 5D gravity in the middle or we have 4D gravity at 
all times. So we require another mechanism to modify early times in a induced gravity 
brane model. 



3.3 Gauss-Bonnet Brane Worlds 

The Gauss-Bonnet term added to the Einstein-Hilbert term, gives the most general action 
in 5D with 2 nd order field equations, as shown by Lovelock [57]. In four dimensions the 
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Gauss-Bonnet term is a topological invariant for compact manifolds without a boundary. 
Therefore Lovelock type gravity in 4D reduces to standard GR. The Gauss-Bonnet terms 
have also been shown to be a low order string correction to the action (see, e.g., 113511361 
W\ El E2 SB SD 02 M M, El HI). In a 5D theory the GB terms are only present in 
the bulk. The bulk action is ll37l: 



Sgrav = 7^2 [ d 5 X^& [R® - 2A £ 
ZK 5 J 



+ a 



R (B)2 _ AR W R (S)ab + R Q^ R (S)<*<A 



(3.11) 



where a is the Gauss-Bonnet coupling. In a classical Gauss-Bonnet theory a can be of 
either sign. It has been shown in Ref. Il58l that for Gauss-Bonnet brane worlds a negative 
value of a gives rise to antigravity or tachyon modes on the brane. However if a bulk 
scalar field is present a negative a can be allowed without these effects occurring. 

The field equations are much more complicated with the Gauss-Bonnet term in the 
action. The Friedmann equation that we obtain is of the form ||39l : 



B , = C ±± C Z -2_K 
8a cr 



where: 



C 



± 



l + -aA 5 + 8a^J + J -±4(p + \)< 



2/3 



(3.13) 



The bulk cosmological constant A 5 has a modified relation to the curvature scale due to 
the GB gravity (55): 



A, = 



6 12a 



(3.14) 



The first term is the standard RS result while the second is the Gauss-Bonnet contribution. 
This also gives us a constraint on a as Eq. (13.141) implies: 



1 1 

& ~ lot 



1 + 3«A 5 



(3.15) 
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We have only one root here as this is the only one that gives the correct RS limit. The 
plus root could be written as: 



4a / 4 

-p=l + Jl + -aA 5 . (3.16) 



The RS limit (a — > 0) is not consistent with this solution. 
From Eq. (|3.15l) we see that: 



a<—. (3.17) 



As J 1 + |aA 5 > 0, this also ensures that A 5 < 0. 

If we consider the case C = = A 5 and p ^> A for a flat brane (K = 0) we see, from 
Eq. (13.121) . that H ~ p 1 / 3 . The scale factor behaves as: 

a <xt 1/{1+w) . (3.18) 

So the Gauss-Bonnet bulk gravity effects the early universe. At late times C/a 4 , p — > 
with A = = A 5 we can show that we end with H = so again we end with a Minkowski 
brane. 

This GB Friedmann equation has no 4D limit: 

GB high energy: H > a~ 1/2 H 2 oc p 2/3 , (3.19) 
GB low energy: H < a~ 1/2 H 2 oc p 2 . (3.20) 

The Friedmann equations for pure DGP and pure GB models with a Minkowski bulk 
are compared in Fig. 13.21 The GB Friedmann equation is a cubic, so has three roots. We 
can show that for p > there is only one real root. This is shown in Fig. 13.21 When 
p = there is the repeated root H = 0. 



3.4 Gauss-Bonnet-Induced Gravity (GBIG) Branes 

In this section and the next chapter we shall consider what happens when we combine 
the induced gravity terms of the DGP model with a Gauss-Bonnet bulk. We will start by 
presenting the model's governing equations in full generalisation. We shall then consider 
the modifications in relation to the pure DGP model. This means we shall only consider 
models with zero brane tension in a Minkowski bulk. In the next chapter we shall extend 
this to the non-zero brane tension in Anti-de Sitter bulk cases. We will see that even in 
the simplest case we get some striking new effects such as a finite density big bang. Some 
of the work in this section was first presented in [60J . 
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Figure 3.2: DGP and GB solutions of the Friedmann equation (H vs p) for a Minkowski 
bulk. (Brane proper time t flows from right to left, with t = oo at p = 0.) 



3.4.1 Field Equations 

The generalised gravitational action contains the Gauss-Bonnet (GB) term in the bulk 
and the induced gravity (IG) term on the brane: 



S, 



grav 



24 



d 5 X y/^ \r(V _ 2A 5 + a \R^ 2 - AR { 2R [ " )ab + R^R^ 



abed 



+ 



d 4 Xy/—g 



5 J brane 



R-^\ 
r 



(3.21) 



where a (> 0) is the GB coupling constant and r (> 0)0 is the IG cross-over scale and A 
is the brane tension. As in the previous models we assume Z 2 symmetry about the brane. 
The standard energy conservation equation holds on the brane, 



p + 3H(1 + w)p = , w = p/p . 



(3.22) 



The modified Friedmann equation was found in the most general case (where the bulk 
contains a black hole and a cosmological constant, and the brane has tension) in Ref. Il62l : 



K 

H 2 - $ + 4 



K nl 
a 2 6 



rH 2 + r^-^(p + X) 



(3.23) 



Note: the cross-over scale used is the one in Ref. [321 and is half that used in Refs. 
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where K is the brane curvature and $ is determined by: 

, A 5 C 

$ + 2a$ 2 = + — , (3.24) 
6 a 4 

where C is the bulk black hole mass. In the correct limits this Friedmann equation reduces 
to forms equivalent to those previously seen. 

In the rest of this chapter we shall only consider a Minkowski bulk (A 5 = 0) without 
a bulk black hole (C = 0), so that $ is a solution of: 

$ + 2«$ 2 = 0. (3.25) 

Equation (|3.25l) has solutions $ = 0, — l/2o;, but here we only consider $ = 0, since 
the second solution has no IG limit and thus does not include the DGP model. We shall 
consider $ = — l/2a in Chapter |U In this bulk we shall consider a spatially flat brane 
(K = 0) without tension (A = 0). Therefore the modified Friedmann equation is given 
by: 

3.4.2 DGP brane with GB bulk gravity: combining UV and IR mod- 
ifications of GR 

By defining the dimensionless variables: 

7 = — , h = Hr, fi = — r = -, (3.27) 
or z d r 

the GBIG Friedmann equation becomes: 

( 1 h 2 + l) 2 h 2 = (h 2 - f ,)\ (3.28) 

while the conservation equation becomes: 

fi' + 3h{l + w)fi = 0, (3.29) 

where a dash denotes d/dr, and h = a' /a. In defining the variable 7 as above, we have 
both the IG and GB contributions within one parameter. This reduces the dimensionality 
of the solution state space. 



(3.26) 
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Figure 3.3: Solutions of the GBIG Friedmann equation (h(fi)) with positive 7 (7 = 0.1). 
On the left is the DGP (+) model and its Gauss-Bonnet corrections, GBIG 1 and GBIG 
2. On the right is the DGP (-) model and its GB correction, GBIG 3. The curves are 
independent of the equation of state w. Brane proper time r flows from right to left, with 
r = 00 at \i = 0. 

Combining Eqs. (|3.28l) and (13.291 ), we find the modified Raychaudhuri equation: 

h > = Mi+w)(h>-») 3 

(7/i 2 + l)(3 7 /i 2 + l)-2(/i 2 - yU )' ^ J 

The acceleration a" / a = h! + h 2 is then given by: 

a" _ h 2 (^h 2 + l)(3 7 /i 2 + 1) - {h 2 - p)[2h 2 - 3(1 + w)n\ 



(7/i 2 + l)(37/i 2 + l)-2(/i 2 -/i) 



(3.31) 



These shall be of use later. 

The GB correction, via a non-zero value of 7, introduces significant complexity to the 
Friedmann equation, which becomes cubic in h 2 , as opposed to the quadratic DGP (±) 
case, 7 = 0, for which: 

h 2 = i{l + 2//± Vl + 4/i} , (3.32) 

as seen in Eq. (13.101) . This additional complexity has a dramatic effect on the dynamics 
of the DGP (+) model, as shown in Fig. 13.31 The contribution of GB gravity at early times 
removes the infinite density big bang, and the universe starts at finite maximum density 
and finite pressure (but, as we show below, with infinite curvature). Furthermore, there 
are two such solutions, each with late-time self-acceleration, marked GBIG 1 and 2 on 
the plots. Since GBIG 2 is accelerating throughout its evolution (actually super-inflating, 
h' > 0), the physically relevant self-accelerating solution is GBIG 1. 
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Figure 3.4: Solutions of the GBIG Friedmann equation (h(/j)) with negative 7 (7 = 
—0.1). On the left is the DGP (-) model and its Gauss-Bonnet corrections, GBIG 1 and 
GBIG 2. On the right is the DGP (+) model and its GB correction, GBIG 3. The curves 
are independent of the equation of state w. Brane proper time r flows from right to left, 
with r = 00 at // = 0. 

In Ref. [58J it is shown that a negative value of a leads to antigravity or tachyon 
modes on the brane. A negative value of a can give regular solutions if a bulk scalar field 
is present. If we allow 7 < the effect of the GB terms is to produce a finite density 
big bang on the DGP (-) branch, Fig. 13.41 The DGP (-) branch is wrapped back on 
itself instead of the DGP (+) branch as in the 7 > case. The DGP (+) branch is now 
modified in a similar manner as the DGP (-) was in the previous case. The consequence 
of this is that the finite density big bang model no longer self-accelerates. We could get 
acceleration at late times by including positive brane tension as we shall see later. In the 
rest of this chapter we shall confine ourselves to 7 > 0. 

The cubic in h 2 , Eq. (13.281) . has three real roots when < 7 < 1/4 (see below). Two 
of these roots correspond to GBIG 1 and GBIG 2, which are modifications of the DGP 
(+) model. The third root GBIG 3 is a modification of the DGP (-) model, as illustrated 
in Fig. l3.3[ Note that the curves in these figures are independent of the equation of state w 
of the matter content of the universe - w will determine the time evolution of the universe 
along the curves, via the conservation equation (13.291) . 

The plots show that GBIG 3 starts with a hot big bang, p = 00, in common with the 
DGP (±) and GB models in Fig.[H By contrast, GBIG 1 and GBIG 2 have no big bang, 
since the density is bounded above: 



(3.33) 
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where y,\ (which is positive only for 7 < 1/4), is found below, in Eq. (I3.48I ). 

The finite-density beginning was pointed out in Ref. [62], where the cubic for the gen- 
eral case (i.e., with brane tension, bulk cosmological constant and bulk black hole) was 
qualitatively analysed. The analysis shows that one solution, GBIG 3, is not bounded, 
which was not noticed in Ref. [62J. The numerical plots of the Friedmann equation in 
Fig. l3.3l are crucial to a proper understanding of the algebraic analysis of the cubic roots. 

A detailed analysis ll62l l63l of the cubic equation (|3.28l) confirms the numerical re- 
sults, and shows that (for /i > 0): 




1 

4 
1 



3 real roots, GBIG 1-3, 



(3.34) 




1 real root, GBIG 3. 



(3.35) 



The real roots are given as follows: 



• For < 7 < 1/4: the roots GBIG 1 



2 are 




(3.36) 



where n = 4 for GBIG I, n = 2 for GBIG 2, and the root GBIG 3 is 




(3.37) 



• For 7 > 1/4: the only real root GBIG 3 is 




1 - 2 7 



+ S+ + 3 



(3.38) 



3 
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In the above, S±, Q, R, 6 are defined by 

1/3 



s. 



± 



R ± ^R 2 + Q 3 



(3.39) 



g = ^(l + 2 /U )-i(2 7 -l) 2 , (3.40) 
R = + y (27 - 1)(1 + 2/i) - 2^7^, (3.41) 



cos 30 = R/^-Q 3 . (3.42) 



The GBIG 1 and GBIG 2 solutions agree with those in Ref. [62J (where the roots are 
given in the fully general case, with brane tension and a bulk black hole and cosmological 
constant). 

The explicit form of the solutions can be used to confirm the features in Figs. l3.3l and 
13.41 Equations (I3.37I) - (I3.42I) show that GBIG 3 starts with a big bang, h, fi — > oo, with 
h 2 ~ /i 2 / 3 near the big bang. This is the same as the high-energy behaviour of the pure 
GB model, as shown by Eq. (13.191 ) - the GB effect dominates at high energies in GBIG 
3. This is not the case for GBIG 1-2, where the high energy behaviour is completely 
different from the pure GB model (and from the DGP(+) model). 

The maximum density feature of GBIG 1-2 is more easily confirmed by analysing 
the turning points of \x as a function of h 2 . The GBIG Friedmann equation (13.281) gives 



dfj, 2(/i 2 -/i)-(7/i 2 + l)(3 7 /i 2 + l) 



d(h 2 ) 2(h 2 - (i) 

Substituting dfx/d(h 2 ) = into Eq. (13.281) . we find that 



(3.43) 



K - (3.44) 

* - 2 37)3/2 • 

The second equation shows that positive maximum density only arises for the upper sign 
and with 7 < 1/4, in agreement with the cubic analysis. 

oo> /!;>(] 0<7<^. (3.46) 
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Thus the initial Hubble rate and density for GBIG 1-2 are 



1 + Jl - 37 

hi = \ (3.47) 

37 



2 - 97 + 2 (1 - 37) 3/2 

~27f 



to = — ■ ^ 3 - 48) 



If j = 0, then GBIG 1-2 reduce to DGP(+), and hi = ft = oo. Note that 

h i >2^H i >-. (3.49) 

r 

The case 7 = 1/4, /x. = corresponds to a vacuum brane with de Sitter expansion, with 
h = hi = 2, generalizing the DGP(+) vacuum de Sitter solution 1(341 . 
The late-time asymptotic value of the expansion rate, as ji — > 0, is 



1 ± VI - 47 

/*oo = ^ (3.50) 

2 7 

where the minus sign corresponds to GBIG 1 and the plus sign to GBIG 2. In the limit 
7^0, GBIG 1 recovers the DGP(+) case, h^ = 1, while for GBIG 2, h^ — ► 00; the 
parabolic GBIG 1-2 curve in Fig. [33] "unwraps" and transforms into the DGP(+) curve. 
Equations (|3.34l) and (13.501) show that 

l<h OQ <2 for GBIG 1, (3.51) 

while: 

2 < h^ < 00 for GBIG 2. (3.52) 
The behaviour of the key GBIG 1-2 parameters is illustrated in Figs. l3.5l and (3T6l 



3.4.3 Cosmological Dynamics 

The GBIG 1 model, which is the physically relevant generalisation of the DGP(+) model, 
exists if Eq. (|3.34l) holds. By Eq. (13.271) . this means that the GB length scale L gh = y/a 
must be below a maximum threshold determined by the IG cross-over scale: 

1 l/3" 
7 < 4 L gh = < -y -r. (3.53) 

If the GB term is taken as the correction term in certain string theories, then L g b ~ -^string, 
while r ~ H X , so that this bound is easily satisfied. 
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Figure 3.5: The dependence in GBIG 1-2 of the initial expansion rate and density on 7. 
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Figure 3.6: The GBIG 1-2 late-time asymptotic expansion rate as a function of 7. 
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When Eq. (13.531 ) holds, the universe starts with a maximum density p\ and maximum 
Hubble rate Hi, and evolves to an asymptotic vacuum de Sitter state: 

r 3 f 12a / 8a^ 3/2 ' 



o<p < ^i^pv-^n 1 -^) >> (3 - 54) 



H^kH < H i = ^ l + A /l-4|. (3.55) 



■ia \ V r 2 

At any epoch t , the proper time back to the beginning is: 

f a ° da 

to-ti = r — . (3.56) 
J a . ah 

Since a and h are nonzero on the interval of integration, the time back to the beginning 
is finite. 

The current Hubble rate can be approximated by the final de Sitter Hubble rate, H ~ 
Hoo, so that by Eq. (|3.51l) . the cross-over scale obeys: 

Ho 1 < r < 2H, 1 . (3.57) 

In the DGP(+) limit, r ~ 1.2/Jq" 1 [|64l . The effect of GB gravity is to allow for increased 
r but not beyond r ~ ZHq 1 . 

However, there is a UV-IR "bootstrap" operating to severely limit the GB effect at 
late times. The key point is that appreciable late-time GB effects require an increase in 7 
(see Fig. 13.61) . whereas the primordial Hubble rate H x is suppressed by an increase in 7 - 
as shown in Fig. 13.51 Equations (13.471) and (13.501) imply that: 

H^Hq 7 < - . (3.58) 

Thus the GBIG 1 model does not alleviate the DGP(+) fine-tuning problem of a very 
large cross-over scale, r ~ Hq 1 ~ (10~ 33 eV) _1 . 
The GBIG 1 Friedman equation (13.361) gives: 



3r 2 1 16a //16a \ 2 64a 2 . o92N L tt 

H =64a^ r-3^ + 2 \ U^- 1 -tpr(* + ***P)«» 



e( P ) + T 



(3.59) 
where: 



cos S9(p) 



2048a 4 /i 2 + 96aV(l + 2p) (f| - l) - 2 (ff - l) : 



3r s 



(#-l) 2 - ^(1-2a*)1 3 2 



(3.60) 
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A more convenient form of the Friedmann equation follows from solving Eq. (13.281 ) for 

fj, = h 2 - h(-fh 2 + 1) , /ioo < h < h . (3.61) 
By expanding to first order in h 2 — h 2 ^, we find that at late times: 

h 2 = hl + 2(-\-)v + 0(v 2 ). (3.62) 



Taking the DGP(+) limit — > 1, and comparing with Eqs. (|3.4t and (|3.32l) . we find that 
the effective Newton constant in GBIG 1 is: 

G = ( -\-) ^ , (3.63) 

where G 5 = nl/8n is the fundamental, 5D gravitational constant. In the DGP (+) case, 
G = G B /r. 

Equation (13.631) gives a relation for the fundamental Planck scale M 5 : 



hr^r — > (3-64) 



rif \ M p 2 
' \2-rH J r ' 

where M p is the effective 4D Planck scale, and we used ~ H . As r — > (its 
upper limit), so M 5 increases. This is very different from the DGP(+) case, where: 

Ml = M 2 /r, (3.65) 

so that M 5 is constrained to be very low, M 5 < 100 Me V. In principle, GB gravity al- 
lows us to solve the problem of a very low fundamental Planck scale in DGP(+)- but in 
practice, the UV-IR bootstrap, Eq. (13.581) . means that 7 ~ so that M 5 is effectively the 
same as in the DGP(+) case. Thus the GB modification of the DGP(+) does not change 
the fine tuning of the cross-over scale r, nor the consequent low value of M 5 . This is 
because 7 is forced to be close to if we want self-acceleration to replace dark energy, as 
in the DGP(+) case. However, no matter how small 7 is, a nonzero 7 dramatically alters 
the early universe, by removing the infinite-density big bang. 

What is the nature of the beginning of the universe in GBIG 1? We can use Eq. (|3.43l) 
in Eq. (13.301) . for matter with w > —1, to analyze the initial state, dfx/d(h 2 ) — > 0+. We 
find that h[ = —00, i.e., infinite deceleration: 

a" = -00 . (3.66) 
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(For GBIG 2, with dp/d(h 2 ) — > 0—, we have h[ = +oo.) The initial state has no 
big bang, but it has infinite deceleration, and thus infinite Ricci curvature. The brane 
universe is born in a "quiescent" singularity. Although similar singularities may be found 
in induced gravity models 11651 . they arise from the special extra effect of a bulk black 
hole or a negative brane tension which is very different from the gravitational GB effect 
that is operating in the GBIG singularity. A key point is that neither the DGP(+) model 
nor the GB model avoid the big bang, as shown in Eqs. (13.51) and (13.191) . But together, 
the IG and GB effects combine in a "nonlinear" way to produce entirely new behaviour. 
If we switch off either of these effects, the big bang reappears. 

This singularity is reminiscent of the "sudden" (future) singularities in General Rel- 
ativity ll66ll - but unlike those singularities, the GBIG 1-2 singularity has finite pres- 
sure. The initial curvature singularity signals a breakdown of the brane spacetime. The 
(Minkowski) bulk remains regular, but the imbedding of the brane becomes singular. 
Higher-order quantum-gravity effects will be needed to cure this singularity. 

By performing an expansion near the initial state, using Eq. (13.611) . we find that the 
primordial Hubble rate in GBIG 1 , after the infinite deceleration at the birth of the uni- 
verse, is given by: 



3Hf - nip 



2„ lV2 



8a (1 + AaH 2 



( A - pf/\ (3.67) 



This is independent of the equation of state w so that the universe decelerates for a fi- 
nite time after its infinite-deceleration birth, regardless of the matter content. If there 
is primordial inflation in the GBIG 1 universe, then the acceleration a" will become 
positive. For a realistic model (satisfying nucleosynthesis and other constraints), a" 
must subsequently become negative, so that the universe decelerates during radiation- 
and early matter-domination. Finally, a" will become positive again as the late universe 
self-accelerates due to the IG effects. 

We can simplify the expression (13.311) for the acceleration in GBIG 1 via Eq. (13.611) : 

_ x(37x + 1 - 2y/x) + 3(1 + w)[xy/x - xj^x + 1)] ^ 
37a: + 1 — 2a/x 

where / = a" /a, x = h 2 . For a given w(x), we can plot f(x). We show an example in 
Fig. 13.71 of a simple model, with primordial inflation followed by radiation domination, 
followed by late-time self-acceleration. We have used the effective equation of state: 

n(hf - /4) + /4 < x < hf , 
w = < (3.69) 

hlo <x<n{h 2 i -hl )+h 2 00 . 
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Figure 3.7: The acceleration / = a" /a vs x = h 2 , for a GBIG 1 cosmology with inflation, 
followed by radiation domination, followed by late-time self- acceleration. Brane proper 
time flows from right to left. Here 7 = 1/12, and n = 0.8 in Eq. (13.691) . 

Here < n < 1 is a parameter determining the time of reheating (with n = cor- 
responding to no inflation and n — 1 to no reheating/ radiation). The values used for 
the equation of state are chosen in order to mimic slow-roll inflation in the first epoch 
(w = —0.9) and radiation in the second (w = 1/3). 

3.4.4 Nucleosynthesis and the age of the universe 

We have noted the UV-IR bootstrap, which enforces a very small value on 7 if the uni- 
verse is to be old enough. In order to make this more quantitative we impose constraints 
from nucleosynthesis and the age of the universe. 

As already stated the DGP cross over-scale must be of the order of the Hubble scale 
or larger, therefore: 



where H is the current Hubble rate. Taking nucleosynthesis to have occurred at an 
energy scale of lMeV and the observed universe to be at a scale of 10 _33 eV we have the 
constraint that: 




(3.70) 



(3.71) 
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We can also put a constraint on 7 from requiring H\ to have occurred at an energy scale 
greater than ITeV. This is in order to put it outside the range of collider experiments. 
Using the solution for h\ in Eq. (13.471) we can write 7 as: 



2/ii - 1 
~3hl 



7 = < 3 - 72 ) 



where: 



Therefore 7 is constrained by: 



h = Hp > 10 45 . (3.73) 



7 < 10- 45 . (3.74) 
Taking r ~ H^ 1 ~ (1CT 33 eV) -1 the above constraint gives: 



a < lO^eV - ^. (3.75) 
The string energy scale (a -1 / 2 ) is then constrained by: 

a -i/2 > io- 10 eV. (3.76) 
This is compatible with the constraint from proton decay which gives lloTl : 

a -l/2 > 1Q 17 e y_ (3 y 7) 

The GBIG Friedmann equation, in the case A = = $, can be written as: 

- 2 N 2 H 2 I 80: \ 2 



where the left hand side is exactly zero in the GR limit. To see if the GBIG model 
approximates GR in the era of nucleosynthesis we write the above Friedmann equation 
in a dimensionless form and apply the constraints above. The dimensionless Friedmann 
equation, Eq. (13.781) . can be written as: 

(i-#) 2 ^( 1+ ^ V ) 2 ' <" 9) 

Expanding out the right hand side and using the constraints above we find: 

1-^) 2 <10" 12 , (3.80) 

during nucleosynthesis. Thus the model will safely meet nucleosynthesis constraints if 
Eq. (I3J3T) holds. 



CHAPTER 3. THE DGP, GB AND GBIG BRANE WORLD COSMOLOGIES 53 

In order to calculate the age of the universe we define a new set of dimensionless 
density parameters: 



16 1 $ K 2 K 2 

= -jaHl fi rc = to* = On = ^P, a = g^A. (3.81) 

In terms of the f2 parameters the (general) GBIG Friedmann equation, Eq. (13.231) . can be 
written as: 



a c {2 + n a (E 2 + a)} 2 (e 2 - 2a) = [e 2 - n m (i + z f - a] 2 , (3.82) 

where E = H/H . 

In order to reduce the number of free parameters in the model we consider the con- 
straint on a at the present time (z — 0, H — H ). This gives: 



0* ± = Ti — 7~FT\ K ± V = = - 2 \. (3.83) 



i ] (i-n m -n x 

With non-zero $ and A, E\ is given by: 



E 2 = {2 - 6fi„a c + 9^a c ^ ± 2^1 - 6nMr c - 9Q 2 a Q T n^ . (3.84) 

By considering dz/dE = in the Friedmann equation we get: 



^ _ ^ a c (4aa + _ 4 ) _ 2^ A + 2 (i - 4a c a) E ; 2 - 3a c *W | 1/3 i 



(3.85) 

We will only consider the case of $ = = A, for which E\ and z- x reduce to: 



E 2 = 

1 9na„,, 



^- {2 - 6aa c ± 2v/i-6aa c } , 0.86) 



and: 

f2(i- 4a a) e 2 - 3a a^ 4 - 4a 1 1/3 

Zi = I rc ^ - I - 1, (3.87) 

with: 

a ± = ± (1 ~g_ m) - 2. (3.88) 

V^r c 
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Figure 3.8: The left plots use the plus branch of Eq. (|3.83l) . the right plots use the negative 
branch. Values used are f2 rc = 0.1, Cl m = 1/4, = and £l\ = 0, so that Vl a+ = 0.372 
and £l a _ = —4.372. The top plots show the Hubble rates against redshift. The bottom 
two plots are the equivalent plots in h, fj,. 
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The DGP limit, fi a — > 0, gives E{ = oo due to the presence of fi Q in the denominator, 
and thus of Zi — oo. 

We see, from Eq. (13.881) . that when tt a± = 0: 



tt m = l±2^Q Tc . (3.89) 
Therefore for tt a± to be positive we require: 



tt m < l-2^/tt Tc , forft a+ , (3.90) 
tt m > 1 + iorfl a _. (3.91) 



Since £l a± > implies a > 0, Eqs. (13.901 ) and (13.911) describe the GBIG 1 model. 
When n a± < 0, we have a < 0, and this corresponds to the non-self-accelerating GB 
modifications of the DGP(-) model, as discussed in section I3A21 
For Q a± to be negative we require: 



tt m > l-2V^r c , {ortt a+ , (3.92) 
O m < 1 + 2^0^, forfi Q _. (3.93) 



We only consider < Q m < 1, therefore is always negative and the solutions do 
not self-accelerate. Q a+ has two regions in the fi m , Q Tc plane. One has solutions that 
self-accelerate (0 Q > 0), one that does not (0 Q < 0). If we included a non-zero 
or f^, self-acceleration would be possible on both branches. Fig. [3J2 illustrates the two 
branches for Q\ = = fl^. 

The initial redshift of the universe must be at least large enough to accommodate 
nucleosynthesis, i.e. z- x > 10 10 . This enforces an extremely small value for Vl a (f2 Q « 
5.5xl0~ 5 ), as illustrated in Fig. [3T9l 

For each solution of z- x there is a region of the Q m , fi lc plane that is significantly 
different from the DGP model. This region, with small f2 m and il Ic , is shown for the 
Zi = 1100 case in Fig. 13.101 This region, which is disallowed by observations, is present 
due to the form of the denominator in Eqs. (|3.87l) and (13.841) . Ei becomes infinite when 
fi fc = and Zi — ► oo for fi m = 0. 

We can work out the age (look-back time) of the universe by evaluating the following 
integral: 

f Zi dz 

to ~ U = H, 1 / — (3.94) 

Jo (l + z)E{z) 

Using H = rSkms^Mpc" 1 (H = 73 ± Skins" 1 Mpc" 1 is the WMAP 3 year data 
result (681) we can get the results in years, Figs. 13.1 ll and 13.121 In Fig. 13.1 II we see a 
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Figure 3.9: The Q m , f2 rc plane for GBIG 1. The central dashed line is for Q a = 0. On 
the left the dashed lines are z\ = 10 and the solid lines are z\ = 2. On the right the solid 
lines are z\ = 1100, the decoupling redshift. The lines above f2 Q = have Q a < 0, while 
those below f2 Q = correspond to Vl a > 0. The magnitude of f2 Q increases as you move 
away from the Vl a = line. 
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Figure 3.11: The age of the universe for Q+ . 

sharp transition along the Q a = ridge. This is because the 0^ > and £l a < branches 
are disjoint. In the region where fi Q > (left hand side) Z\ is for the self-accelerating 
branch. On the right hand side of the plot Q a < 0, so z\ is for the non-self-accelerating 
branch. From Fig. 13.121 we see that in general, solutions that lie in the fi Q < region 
will give us larger ages. If we considered solutions with a matter density less than that 
observed (fi m « 0.3) it is possible to get a sufficiently large ages with f2 rc measurably 
different from that of the DGP (0 Q = 0) model. If we want to have self-acceleration 
and the finite density big bang (Q a > 0) and have enough time for the evolution of the 
universe we are restricted to solutions that lie along (or at least extremely close to) the 
DGP limit. This means that observationally the GBIG model will be indistinguishable 
from the DGP model even if the early universe is dramatically different. 

In Figs. 13.131 and 13.141 we have the age results for the case. As f2 a _ < for all 
values for fi m < 1 all the results are for the non-self-accelerating solution. 

3.5 Conclusions 

We have seen that by including the Gauss-Bonnet bulk term with the Induced gravity 
brane term in the gravitational action, we get some striking new features. The ultra-violet 
correction from the Gauss-Bonnet term in combination with the infra-red correction from 
the induced gravity term gives rise to a solution that starts with a finite density and ends 
with self-acceleration. We saw that the infra-red induced gravity term on its own can 
give a self-accelerating solution which is still present in the GBIG model. The GB term 
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Figure 3.12: As in Fig. 13. 1 ll but viewed from above. 
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Figure 3.14: As in Fig. I3.13L but viewed from above. 

modifies the early era but but does not produce the finite density big bang on its own. It 
is only the combination that does. There is no threshold for the finite density big bang to 
occur- the contributions from the GB and IG terms only need to be non-zero. 

The density, pressure and temperature are finite but the big bang does have a curvature 
singularity at the beginning. This deserves further investigation. 

We have seen that the UV-IR bootstrap confines 7 to be very small i.e. the GBIG 
model is very close to the DGP(+) model for most of the history of the universe. This 
was also seen in the analysis of the age of the universe. In order for there to sufficient time 
for events such as nucleosynthesis to have taken place 7 must be almost indistinguishable 
from the DGP(+) model (7 = 0). 



Chapter 4 



Generalised cosmologies with Induced 
Gravity and Gauss-Bonnet terms 



4.1 Introduction 

In this chapter we shall extend the work in Chapter [3] to models with non-zero brane 
tension and bulk curvature. We keep the assumptions that the bulk black hole mass is 
zero and the brane is spatially flat. We shall see that this leads to many more possible 
solutions, but not all of them are physically relevant. 



4.2 Field Equations 



The general form of the Friedmann equation in this case is: 



1 + -a I H 2 + ' 
3 V 2 



(H 2 - $) 



rH 2 -^(p+X) 



where $ is a solution to: 



(4.1) 



$ + 2a<£> z 



As 
6 



(4.2) 



Due to the presence of the GB term in the bulk action, the bulk cosmological constant is 
given by: 



As = + 



6 12a 



£ 2 £ 4 ' 

as seen in Chapter[3] Equations (14.21) and (14.31) give us two solutions for $: 



(4.3) 
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$1 is the generalised form of the solution considered in the previous chapter. 

We work with the Friedmann equation in dimensionless form. We shall use the same 
variables as defined in the previous chapter (Eq. 13.271) with the addition of: 



2 2 

rni , r 



o = -f\ X = ^,<P = (4.5) 

Using these variables we can write the two solutions for $ as: 

4 

37 

The bulk cosmological constant gives us an upper bound on the GB coupling constant 
a. Equation (14.31) gives us: 



1 

4a 



I ± \/L + ~aA 5 



(4.7) 



For a RS (a — > 0) limit we take the minus branch. A 5 > would imply £ 2 < and thus 
we must have A 5 < and: 



I 2 

a<—. (4.8) 



In dimensionless form this is given by: 



7 < i (4-9) 

Maintaining a RS limit would also rule out the 2 branch. We would be restricted to the 
solutions lying along the line in the top left quadrant in Fig. 14. II We are interested in the 
whole range of the model so we include the plus branch in Eq. (I4.7I ). We assume A 5 < 
therefore the constraint on a is given by: 

e 4 

a<-, =>- 7 < 7T- < 4 - 10 ) 

2 3x 

If we take x — then we have no bound on 7 (apart from being positive and real). This 
is the case considered in Chapter[3l 

In Fig. 14. II we have the two solutions plotted as functions of x- The two solutions 
with = both live in a Minkowski bulk, all the rest live in an AdS bulk. Note that one 
of these AdS solutions (% = 0, = —4/37) nas A 5 = but $ = —l/2a acting as an 
effective cosmological constant, Chapter [3] We see that for any allowed value of we 
can be on either of the two branches. This means that we need not consider the <fii and 
02 solutions to the Friedmann equation separately. We therefore consider Eq. (|4.10l) in 
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Figure 4.1: The effective bulk cosmological constant as a function of x- 



terms of 0. We define the maximum value of 7, for a particular value of 0, that is allowed 
by the constraint equation: 



7m = - 



The dimensionless Friedmann equation is: 



30' 



(4.11) 



with the conservation equation now given by: 



(4.12) 



/i' + 3A(1 + u))/j = 0, 



(4.13) 



where ' = d/dr and h = a' /a. The Raychaudhuri and acceleration equations are given 
by: 



ti 



3fi{l + w)[h 2 - {fi + a)] 



(7/i 2 + l)(3 7 /i 2 + 1) - 2[h 2 - (fj, + a)] - 07 (1 + |0 7 ) ' 



(4.14) 



and: 



1" h 2 ^ 2 + l)(3 7 /i 2 + 1) - [h 2 -(/! + a)][2h 2 - 3/i(l + w)] - /i 2 07 (l + |07) 



(7/1 2 + l)(3 7 /i 2 + 1) - 2[h 2 + a)} - 07 (1 + |0 7 ) 



(4.15) 
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Figure 4.2: Solutions of the Friedmann equation (h vs fi) with negative brane tension 
(er = —4) in a Minkowski bulk (0 = 0) with 7 = 1/10. The curves are independent of 
the equation of state to. The arrows indicate the direction of proper time on the brane. 

4.3 Friedmann Equation Solutions 

The model we considered in Chapter [3] is the Minkowski bulk limit (x = 0) °f me 0i 
case, which we now see to be equivalent to % — j- m me ^ 2 case - The results we present 
below are in terms of and thus include both 0i and 02. We will consider the effect of 
including brane tension in a Minkowski bulk before we look at the AdS bulk cases. The 
Minkowski bulk case is given by = 0, since Eqs. (|4.3I) and (|4.2I) imply A 5 = 0. 

4.3.1 Minkowski bulk (0 = 0) with brane tension 

In Figs. l4.2l and 14.31 we can see the solutions to the Friedmann equation in a Minkowski 
bulk for both positive and negative brane tensions. By introducing brane tension into 
the model we are effectively adding a brane cosmological constant. These solutions 
are therefore less desirable than the zero brane tension, self-accelerating ones. For both 
negative and positive brane tensions there are three solutions (this is not always true as we 
will show later) denoted GBIG 1-3. There are four points of interest in Figs. l4.2l and l4.3[ 
these are: 

• (/ii, h{) and (/i e , hi): The initial density for GBIG 1 and GBIG 2 (/ij) and the final 
density for GBIG 1 and 3 (/x e ), found by considering d^/d(h 2 ) = 0, are given by: 
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Figure 4.3: Solutions of the Friedmann equation with positive brane tension (a = 2) in a 
Minkowski bulk (0 = 0) with 7 = 1/10. 



where the plus sign is for /i; and the negative sign is for /i e . The Hubble rates for 
these two densities are given by: 

V = i±^35, (4.17, 

37 

where the sign convention is the same as above. The points h he ) have h' = 
|oo|. Therefore cosmologies that evolve to /i e , h e end in a "quiescent" (finite den- 
sity) future singularity. This singularity is of type 2 in the notation of Ref. Il65l . 



(/i/, 0): This is the density at which GBIG 3 "loiters". The point was found by 
considering h = and is given by : 



IH = -a. (4.18) 

In the case considered in Chapter [3] we had a = so fii = 0. We can show that 
for a < GBIG 3 will not collapse but will loiter at a density of /i = —a before 
evolving towards (fi e , h e ), by considering h\. At = —a, h = we get h\ = 
from Eq. (14. 141 ). In a standard expanding or collapsing cosmology h! < at all 
times. The evolution for radiation and dust dominated universes can be seen in 
Fig. 14.41 A dust dominated universe loiters for longer than the radiation dominated 
universe. The length of time a universe will loiter for is also dependent on 7 in a 
non-trivial way. This dependence is a matter for further investigation. In Ref. ||69l 
they consider a loitering braneworld model. An important difference between the 
two models is that in Ref. ||69l they require negative dark radiation, i.e. a naked 
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Figure 4.4: Plots of fj, vs r and h vs r for GBIG 3 in a Minkowski bulk ((f) = 0) with 
negative brane tension (a = —1/4). We see that GBIG 3 loiters around fi = Hi = — a, 
h = before evolving towards a vacuum de Sitter solution (/i e < 0). (Here 7 = 1/10. 
The solid line is w = 1/3. The dotted line is w = 0.) 



singularity in the bulk or a de Sitter bulk. Also in the model the Hubble rate at the 
loitering phase is exactly zero. The time spent at this point is only dependent on w 
and the density of the loitering phase is only dependent on the brane tension. 

(0, /too): This is the asymptotic value of the Hubble rate as \i — > 0. For the value of 
the parameters used in Fig. 14.21 GBIG 2 is the only case with this limit. In general 
any of the models can end in a similar state (Fig. |4.3l ). The different values of 
are given by the solutions to the cubic: 

/t 6 + (27 ~ X V + ^h* - - = (4 19) 

In the case we considered in Chapter [3] (a = 0) the last term in the above cubic is 
zero. Therefore we have = for GBIG 3 while the solutions for GBIG 1 and 
GBIG 2 are given by: 



»- = l±£EE., (4.20) 

27 

where the minus sign corresponds to GBIG 1 and the plus sign to GBIG 2. This is 
the only case where we can write simple analytic solutions as in all other cases we 
have to solve the cubic. 



The two Hubble rates, h\ and h e , are independent of the brane tension which simply 
shifts the [i = axis. Therefore /i c only corresponds to a physical (positive) energy 
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Figure 4.5: Solutions of the Friedmann equation (h vs fi) in a Minkowski bulk (0 = 0) 
with 7 = 1/7. The vertical lines represent the ji = axis for the labeled brane tensions. 



density when a < a c < 0, see Eq. (I4.21I ). The effect of the brane tension on the fi = 
axis is illustrated in Fig. 14.51 

We showed in Chapter |3] that for <\> = = o we require 7 < 1 /4 for GBIG 1 and 
GBIG 2 solutions to exist within the positive energy density region (if 7 = 1/4 GBIG 1 
and GBIG 2 reduce to the same vacuum de Sitter universe). If we have some non-zero 
brane tension this constraint is modified. The maximum value of 7, for GBIG 1 and 
GBIG 2 to exist with positive energy density, as a function of a can be seen in Fig. 14.61 
We have defined new quantities, cr c i and 07, for which /i e i 1 = 0. We can see from 
Eq. (14.181) that 07 = 0. cr e j are given in terms of 7 by: 

2-9 7 ±2(l-3 7 )'/' 
277 2 

The maximum value of 7 (7 m ) for GBIG 1 and GBIG 2 to exist (i.e. for Eqs. (14.161) 
and ( 14.171) to have real solutions) is 7 m = 1/3 . Actually for 7 = 1/3, GBIG 2 exists 
but GBIG 1 is lost. This is since hi = h e at 7 = 1/3. The point hi = h c is now a point 
of inflection and Eq. (14.171) is still valid and therefore \h'\ = 00. For values of 7 > 1/3, 
hi >e and yU; e become complex, it is no longer a point of inflection. Therefore \h'\ 7^ 00 
and GBIG 3 can continue its evolution through this point to h^. In Fig. 14.71 we show 
results for h and h! for three values of 7, 7 = 1/3 — 0.01, 1/3, 1/3 + 0.01. In the top 
plot we see that the solution denoted (1) has GBIG 1-3 present as 7 < 1/3, we have two 
real and different values for h- x and h c . This solution in the bottom plot has three parts 
(GBIG 1 has negative values not shown in Fig. 14.71) . GBIG 2 is the dotted solution and 
the solid is GBIG 3. As (1) approaches h lfi \h'\ — > 00. Solution (2) has 7 = 1/3, in 
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Figure 4.6: The (a, 7) plane for solutions in a Minkowski (0 = 0) bulk. The short dotted 
horizontal line is 7 m = 1/3. GBIG 1-2 exist with positive energy density in regions III, 
IV and V. 
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Figure 4.7: The plot on the left is h vs /x for three different solutions, all with = 0. The 
plot on the right shows h! vs /i for the same solutions. (1): 7 = 1/3 — 0.1, a — —3. (2): 
7 = 1/3, a = —2.5. (3): 7 = 1/3 + 0.1, a = —2. Different brane tensions are used 
purely for clarity. 
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the top plot we see that GBIG 1 now no longer exists as hi = h c . In the bottom plot 
GBIG 3 (solid) and GBIG 2 (dotted) solutions both go to h' — oo at the point hi = h e . 
Solution (3) has 7 > 1/3. Now h iie does not represent dfx/dh = 0, therefore h! stays well 
behaved throughout the evolution (bottom plot). GBIG 3 moves seamlessly onto GBIG 2 
making a single solution. As GBIG 2 super-accelerates this new solution will show late 
time phantom-like behaviour ih! > — > w < — 1). 

Solutions that lie along the 07 line in Fig. 14.61 are those considered in Chapter [3] 
Regions I and II in Fig. 14.61 extend up to 7m = 00. Solutions in each region of Fig. 14.61 
have the following properties: 

• I : a c (7 m ) < a < 0, 7 > 7i and a < cr e (7 m ), 7 > 7 m . GBIG 1 and GBIG 2 do 
not exist. GBIG 3 loiters at fii before evolving to a vacuum de Sitter universe. 

• II : a > 0, 7 > 7i. GBIG 1 and GBIG 2 do not exist. GBIG 3 evolves to a 
vacuum de Sitter universe. 

• III : a < a e , 7 < 7 m . GBIG 1 will evolve to (// e , h e ). GBIG 2 evolves to (0, hoc). 
GBIG 3 loiters at fii before evolving towards (/i e , h e ). When a = a c GBIG 1 and 3 
evolve to (0, h c ). When 7 = 7 m , GBIG 1 ceases to exist (hi = h e ). 

• IV : a e < o < 0, 7 < 71. GBIG 1 and GBIG 2 both evolve to vacuum de 
Sitter states. GBIG 3 loiters before evolving to a vacuum de Sitter universe. Each 
vacuum de Sitter state has a different value of h^. When 7 = 7; GBIG 1 and GBIG 
2 live at (0, hi), 

• V : a > 0, 7 < 7;. GBIG 1-3 all end in vacuum de Sitter states. With 7 = 7; 
GBIG 1-2 live at (0, hi). When a = GBIG 3 ends in a Minkowski state. 

In Fig. 14.81 we can see how solutions in each of the regions mentioned above affect 
/loo. The thin-dark line (a = -0.4) lies in III and I in Fig. glO Only GBIG 2 exists 
for < 7 < 1/3. For 7 > 1/3 we've seen that GBIG 3 and GBIG 2 connect to make 
a single solution, which is why the line is continuous through 7 = 1/3. The thin-light 
line (cr = —0.28) has an interesting feature due to the solutions lying on a line that cuts 
through both region III and IV as well as I. Both the thick lines have GBIG 1-3 ending 
at h^. The thick-dark solution loiters before evolving to h^, but this has no effect upon 
the results in Fig. 14.81 

4.3.2 AdS bulk (0 ^ 0) with brane tension 

When (f) = 0, we have A 5 7^ 0, with one exception: the 2 solution with \ — 0, i.e. 
<p2 = —4/37, has A 5 = 0, but the bulk is AdS. For x > 0, A 5 7^ 0. Thus in all cases, 
<p 7^ implies an AdS bulk. 
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Figure 4.8: for solutions in a Minkowski (0 = 0) bulk. The thin-dark line has 
a = —0.4, thin-light line has a = —0.28, thick-dark lines have a = —0.2 and the 
thick-light lines have a = 0.2. 

When we allow the bulk to be warped (0 ^ 0) we open up another possible solution, 
denoted GBIG 4. There is a maximum value of </> for which GBIG 4 can exist as we shall 
see later. The nature of GBIG 3 is also changed. These solutions, for a negative brane 
tension, can be seen in Fig. I4.9I Brane tension affects the solutions in the same manner 
as in 4> = case. The qualitative effect of the warped bulk is to make GBIG 3 collapse 
and to introduce the new bouncing branch GBIG 4. This is due to h = now giving two 
solutions: 



where the plus sign is for the GBIG 3 collapse (/i c ) and the minus for the GBIG 4 bounce 
(/^b)- Effectively the loitering point in the Minkowski bulk is split into the max/min den- 
sities of the bouncing/collapsing cosmologies of GBIG 3-4. The other points in Fig. |4.9l 
are modified by the warped bulk; they are now given by: 



with the plus sign for /i ; and the negative sign for // e . The Hubble rates h i c are given by: 




(4.22) 



4 - I87 + 277 2 ± y/2(2 - 67 - 97 2 0) 3 / 2 



a. 



(4.23) 



h i>e = ^1/4 - 67 + 9 7 2 ± 2 v / 4^ 12 7 - 18 7 2 0, 



(4.24) 



with the same sign convention, is given by a similarly modified equation: 
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Figure 4.9: Solutions of the Friedmann equation (h vs jj,) with negative brane tension 
(a = —3/4) in an AdS bulk (0 = —0.01) with 7 = 1/5. The curves are independent of 
the equation of state w. The arrows indicate the direction of proper time on the brane. 



ht + 



(2 7 - 1), 4 [1 + 2(7-07(1 + 107)] 



r 



hi 



r 



0(i + ^) 2 + ^ 



0. 



(4.25) 



With a warped bulk the constraint from hi (7 < 1/3 in the Minkowski case) is mod- 
ified. In the equation for hi we effectively have two bounds from the two square root 
terms. In the = case only the inner term is of any consequence (giving rise to the 
bound 7 < 1/3). When 7^ there are two bounds which are applicable in different 
regimes. If we consider the bound from the inner square root we get: 



7< 



^TT20- 1 
30 



(4.26) 



This is valid for > —1/2. Considering the outer square root term we get the bound 

2 (l - 2v^0) 



7< 



30 



(4.27) 



This bound becomes negative when > —1/4 which is disallowed as this prevents 
self-acceleration. The bound on 7 is given by the lower of the two constraints when 
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Figure 4.10: In order for h x to be real 7 must lie beneath the solid line (7m). The region 
to the right of the vertical dotted line is where GBIG 4 (h c > 0) is allowed (except where 



— 1/2 < <p < — 1/4 i.e. when in the range where both exists. Therefore the bound on 7 
for h\ to be real is given by (see Fig. 14. 101 ): 



For a particular value of the maximum value of 7 allowed by this constraint is denoted 
7 m . When < -9/4 the constraint in Eq. (14.101) is tighter than that in Eq. (14.271) (7 M < 
7 m ). This means that for < — 9/4, hi is always real for allowed values of 7. 

There is a bound for GBIG 4 to exist, found by considering h e = 0. This bound is 
given by: 



This is a solution to the quadratic obtained from h e = 0, the other root of the quadratic 
does not obey the constraint < —4/37 so it is ignored. The minimum value for GBIG 4 
to exist will be denoted 4>GBiG4Lim- Using the constraints we can split the 7, plane into 
three sections, see Fig. 14.101 The solid line is constructed from the bounds in Eq. (14.281) . 
For hi to be real we must choose values below this line. The area to the right of the 
vertical dotted line (but excluding = 0) allows GBIG 4. Points to the left of this line 
have GBIG 1 collapsing after a minimum energy density (/i b ) is reached. The dotted 
curve on the left comes from initial constraint 7 < —4/30. 



= 0). 




(4.28) 




(4.29) 
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Figure 4.11: The (a, 7) plane for solutions in a AdS (0 = —0.1) bulk. The short dotted 
horizontal line is the maximum value of 7 as obtained from Eq. (14.101) . 



In order to consider the a, 7 plane as we did in the = Minkowski case, we need 
the modified a^ c equations: 



<Ti,e = 1 ' ^ 1 *^>—. (4.30) 



4 - I87 + 277 2 ± y/2(2 - 67 - 97 2 * 
547 

We now have two more formulas for the collapse density of GBIG 3 and the bounce 
density of GBIG 4: 

a c>h = ±^(l + ?P\ , (4.31) 



2 t 

with the plus sign corresponding to the collapse and the minus to the bounce. 

We shall consider the a, 7 plane for three different values of <fi corresponding to three 
distinct regions in Fig. 14. 101 We shall first consider <\> = —0.1. 

Typical example </> = —0.1 

If we take = —0.1, we are in the region where GBIG 4 is allowed. The a, 7 plane can 
be seen in Fig. 14.1 ll The regions I, II and III in Fig. 14.1 1 1 extend up to 7m = 40/3. 
In each region we have the following cosmologies. 

• I : cr c (7 m ) < cr < 0, 7 > 71 and a < cr e (lm), 7 > 7m- GBIG 1-2 do not exist. 
GBIG 3 reaches a minimum energy density (/i c ) and then evolves back to /i = 00. 
GBIG 4 starts at (0, — h^) then bounces at (/ib, 0) before evolving to (0, +h 00 ). 
When a = o"b, GBIG 4 exists as a Minkowski universe. 
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• II : o"b < cr < cr c , 7 > 7i- Only GBIG 3 exists. GBIG 3 reaches a minimum 
energy density (// c ) and then evolves back to \i = oo. When a = a c GBIG 3 ends 
in a Minkowski universe. 

• III : a > (T c , 7 > 7;. Only GBIG 3 exists and evolves to a vacuum de Sitter 
universe. 

• IV : a < a c < 0, 7 < 7 m . GBIG 1 evolves to (ju e , h e ). GBIG 2 evolves to h^. 
GBIG 3 reaches a minimum energy density (/i c ) and then evolves back to fi = 00. 
GBIG 4 starts at (// e , — h e ), bounces at (/i b , 0) before evolving to (ji e , +h c ). When 
o = a c GBIG 1 and 3 evolve to (0, h e ). GBIG 4 starts at (0, —h c ) and bounces 
before evolving to (0, h e ). When 7 = 7 m GBIG 1 ceases to exist (hi = h e ). 

• V : cr c < cr < (Tb, 7 < 7i- GBIG 1 and GBIG 2 both end in vacuum de Sitter 
universes with different values of h^. GBIG 3 reaches a minimum energy density 
(/i c ) and then evolves back to fi = 00. GBIG 4 starts at (0, — h^) then bounces 
at (/z b , 0) before evolving back to (0, +h 00 ). When a = a h GBIG 4 exists as a 
Minkowski universe (0, 0). 

• VI : Ob < cr < cr c , 7 < 7;. GBIG 1 and GBIG 2 both end in vacuum de Sitter 
universes with different values of h^,. GBIG 3 reaches a minimum energy density 
(/j, c ) and then evolves back to fi = 00. GBIG 4 does not exist. When a = a c GBIG 
3 ends in a Minkowski universe. When 7 = 7i GBIG 1 and GBIG 2 live at (0, hi). 

• VII : a > cr c , 7 < 7j. GBIG 1-3 all end in vacuum de Sitter universes with 
different values of h^. GBIG 4 does not exist. When 7 = 7; GBIG 1 and GBIG 2 
live at (0, hi). 

In region I there is a GBIG 4 solution, which is modified in the same way as the 
GBIG 3 solution in the Minkowski bulk. As h- he are no longer valid in region I GBIG 4 
joins onto GBIG 2, see Fig. 14.121 The plot on the right in this figure shows h", in order 
to clearly distinguish this solution from the one in the = — 1 case where GBIG 4 is no 
longer allowed. 

In Fig. 14.131 are the results for h^ with <p = —0.1. The thin-dark line (a = —1/2) 
lies in regions I and IV. Therefore we have only one solution for h^, which corresponds 
to GBIG 2 for 7 < 7 m . For 7 m < 7 < 7m this root now corresponds to the end 
point for GBIG 3 as in the = case. The thin-light lines (a = —0.35) lie in regions 
I and V. There are the two GBIG 1 and GBIG 2 solutions and the GBIG 4 solution 
which converges with the light-thick line at the bottom. The thick-dark line (cr = 0) 
lies in regions II and VI, so only has GBIG 1 and GBIG 2 present. The thick-light lines 
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Figure 4.12: The GBIG 4 solution in region I for <p = —0.1, 7 = 0.4 > 7 m and a = 
— 1/2. The plot on the right is h" vs // for this solution. This shows the differences 
between this solution and the similar case when = — 1, where GBIG 4 is no longer 
allowed. Arrows denote proper time. 

(a = 1 /2) lie in regions III and VII, so has GBIG 1 and GBIG 2 and the GBIG 3 solution 
(the horizontal line). 

Typical example: <p = — 1 

We now consider = — 1 which is in the region where GBIG 4 no longer exists (see 
Fig. 14.141) . as h e = 0. Therefore fi e is not relevant and GBIG 1 bounces at /ib- This 
means that the (a, 4>) plane, Fig. 14. 151 is simpler. The regions in Fig. |4.15| are: 

• / : cr < cr b , 7m < 7 < 7m- GBIG 1 does not exist. GBIG 2 starts at (0, —h^), 
collapses to (/i b , 0) and then expands back to (0, h^). GBIG 3 expands to (/i c , 0) 
and then collapses. When a = GBIG 2 exists as a Minkowski universe. 

• II : cr b < cr < ^c, 7i < 7 < 7m- GBIG 1 and GBIG 2 do not exist. GBIG 3 
expands to fx c and then collapses. When a = a c GBIG 3 evolves to a Minkowski 
universe. 

• III : o- > o- c , 7; < 7 < 7m- GBIG 1 and GBIG 2 do not exist. GBIG 3 evolves 

tO /loo. 

• IV : a < o"b, 7 < 7m- GBIG 1 evolves from /i; to /ib and then expands back to yU;. 
GBIG 2 evolves to h^. GBIG 3 expands to /i c and then collapses. When 7 = 7 m 
GBIG 1 ceases to exist, GBIG 2 expands from (/i b , 0). When a = a h GBIG 1 
evolves to a Minkowski universe. 
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Figure 4.13: for solutions in a AdS (0 = —0.1) bulk. The thin-dark line has a = 
— 1/2, thin-light lines have a = —0.35, thick-dark line has a = and the thick-light lines 
have u = 1/2. 




Figure 4.14: Solutions of the Friedmann equation (h vs /j) with negative brane tension 
(a = -2) in an AdS bulk (0 = -1) with 7 = 1/10. For this value of GBIG 4 no longer 
exists and GBIG 1 can collapse. The curves are independent of the equation of state w. 
The arrows indicate the direction of proper time on the brane. 
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Figure 4.15: The (er, 7) plane for solutions in a AdS (0 = —1) bulk. The short dotted 
horizontal line is the maximum value of 7 as obtained from Eq. (14.261) . The top horizontal 
line is from the initial bound in Eq. (14.101) . 

• V : a h < a < a c , 7 < 71. GBIG 1 and GBIG 2 evolve to h^. GBIG 3 expands 
to /i c and then collapses. When 7 = 7^ GBIG 1 and GBIG 2 exist as the same de 
Sitter universe with (0, h^). When a = a c , GBIG 3 ends as a Minkowski universe. 



• VI : cr > cr c , 7 < 7 ; . GBIG 1-3 evolve to h^. When 7 = 7;, GBIG 1-2 live at 

(0,/ioo). 



In region I we again have a combined solution as GBIG 1 has vanished. As GBIG 4 
is not allowed (h c = 0) when we take 7 > 7 m , which causes hi = 0, GBIG 2 matches up 
with its negative counterpart. We see the bouncing GBIG 2 solution in Fig. 14.161 Note 
the nature of h" is very different to that of the GBIG 2, 4 bounce in the <p = —0.1 case. 
In Fig. |4.17| we present results for for = — 1. 

Typical example: <p = — 5 

Here we shall present results for = — 5 for completeness. This value of lives in 
the region of Fig. 14.101 where hi is always real. This means that the a, 7 plane is much 
simpler, Fig. 14.181 

The regions in Fig. |4.18| are: 

• / : a < at,, 7 < 7m- GBIG 1 evolves from fii to fib and back. GBIG 2 evolves to 
h^. GBIG 3 expands to /i c and then collapses. When a = a^, GBIG 1 ends in a 
Minkowski universe. 
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Figure 4.16: The GBIG 4 solution in region I for <p = —1, 7 = 0.8 > 7 m and a — — 1. 
The plot on the right is h" vs ji for this solution. This shows the differences between this 
solution and the case when = — 0.1. Arrows denote proper time. 
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Figure 4.17: for solutions in a AdS (0 = —1) bulk. The thin-dark line have a = —1.2, 
thin-light lines have a = 0.8, thick-dark line have a = and the thick- light lines have 
a = 1.2. 
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Figure 4.18: The (cr, 7) plane for solutions in a AdS ((f) = —5) bulk. The top horizontal 
line is the bound from the initial bound in Eq. (14.101) . 

• II : (T; < < (T c , 7i < 7 < 7m- GBIG 1 and GBIG 2 do not exist. GBIG 3 
expands to /z c and then collapses. When a = a c GBIG 3 evolves to a Minkowski 



• III : cr > (T c , 7i < 7 < 7m- GBIG 1 and GBIG 2 do not exist. GBIG 3 evolves 



• IV : Ob < cr < <7i and <Xb < cr < a c ; 7 < 7m and 7 < 7;. GBIG 1 and GBIG 
2 do not exist. GBIG 3 expands to /i c and then collapses. When 7 = 7i (and for 
a = ad GBIG 1 and GBIG 2 both live at (0, h^). When a = a c GBIG 3 evolves 
to a Minkowski universe. 

• V : a > cr c , 7 < 7;. GBIG 1-3 evolve to h^. When 7 = 7;, GBIG 1-2 both live 

at (0, /loo). 

As we decrease the value of <p, region II in Fig. 14.181 shrinks (the point cr;(7 M ) be- 
comes increasingly positive). For (f> < —64/9 region II no longer exists. The nature of 
the solutions in the other regions are unaffected. 

There are no bouncing solutions in this case (in fact for any case with <p < —9/4) as 
GBIG 4 is unallowed and 7 M < 7 m (which rules out the combined solutions). In Fig. |4.19l 
we show the h^, results for <p = — 5. 
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4.4 Conclusions 

In this chapter we have looked at the general GBIG model. We have seen that there is 
a range of possible dynamics that can be achieved depending on the parameters in the 
model. The main model of interest is still GBIG 1 as this is the one that starts in a finite 
density "quiescent" singularity and evolves to a de Sitter universe. If we warp the bulk 
enough, GBIG 1 can be allowed to collapse back to its initial density, provided there 
is sufficient negative brane tension. If the bulk is warped enough to allow GBIG 1 to 
collapse this means that the bouncing cosmology GBIG 4 can no longer exist. The exact 
nature of the late time dynamics can be changed by including a non-zero brane tension. A 
negative brane tension can reduce the Hubble rate at late time and a positive tension will 
increase it. GBIG 1 can end in a future "quiescent" singularity with a non-zero density 
if we have an appropriate (negative) brane tension present. As the solution of interest is 
GBIG 1 and we want it to provide the late time acceleration that we are experiencing it 
needs to end as a vacuum de Sitter universe. Therefore if there is some brane tension it 
must take values o\ > a > a c . 

GBIG 2 also starts in a "quiescent" singularity but it super-accelerates so it is un- 
physical and of little interest. 

GBIG 3 still starts in an infinite density big-bang but has a number of possible late 
time dynamics. In a Minkowski bulk GBIG 3 can either evolve to a Minkowski state, 
as shown in Chapter [3l a vacuum de Sitter state or even loiter around ji = —a before 
ending in a vacuum de Sitter state or a "quiescent" singularity. This loitering cosmology 
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is different from that in Ref. (691, as we do not require a naked bulk singularity or a de 
Sitter bulk. If we warp the bulk then GBIG 3 will generally collapse, unless there is 
sufficient (positive) brane tension to allow the solution to end in a vacuum de Sitter state. 

GBIG 4 can only ever exist in a mildly warped bulk with negative brane tension. So 
it can be said that GBIG 4 is unphysical due to the requirement that a < 0. 

There are a number of bouncing cosmologies within this set-up, in addition to the 
GBIG 4. There are also the solutions where GBIG 4 and 2 match up (j m < 7 < 7m 
with 0GBiG4Lim < < 0) and the GBIG 2 bouncing cosmologies (7 m < 7 < 7m 
with — 9 < < 0GBiG4Lim)- Each of these have different dynamics so would produce 
different evolutionary histories. The solutions that spend time on the GBIG 2 branch will 
experience phantom-like behaviour during this period. 



Chapter 5 
Conclusions 



In this thesis we have looked at some brane world cosmological models. The Randall- 
Sundrum is an interesting toy model of the universe in which we are able to investigate 
some of the phenomenological properties of string theory ideas. The idea that everything 
apart from gravity is confined to a four dimensional hypersurface leads to some striking 
new features. 

In chapter|2]we looked at the Kaluza-Klein modes of the graviton in Randall- Sundrum 
models. We considered the nature of these modes for both Minkowski and de Sitter 
branes in the Randall-Sundrum one and two-brane models. For two Minkowski branes 
there is a zero mode and then a series of massive modes. The spacings of these modes 
are governed by Bessel functions, dependent on the AdS length scale of the bulk and 
the inter-brane distance. When we send the second brane out to infinity we obtain a 
continuum of massive Kaluza-Klein modes. In the investigation of two de Sitter branes 
we saw that there is a zero mode and a mass gap, which is a function of the Hubble rate on 
the brane, to the first massive mode. We derived a new expression for the mode spacing 
as a function of associated Legendre functions, dependent on the AdS length scale and 
the inter-brane distance but also on the Hubble rate on the brane. Therefore they differ 
in the low energy and high energy limits. Taking the one brane limit we again find a 
continuum of massive modes. 

In chapter[2]we reviewed the 5D bulk equations and how these are projected onto the 
brane in Randall-Sundrum models. It is the projected equations that are used when we 
want to understand the cosmological dynamics in these models. 

In chapter |3] we briefly reviewed the DGP and GB models before combining the two 
into GBIG model. The DGP is a very interesting brane world model that much work has 
been done on, because of the late time acceleration in the DGP(+) branch. This is a very 
interesting result as the observed acceleration of the universe is usually interpreted via 
the introduction of a "dark energy" field. The DGP model explains this phenomenon via 
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Figure 5.1: The three GBIG models with zero brane tension in a Minkowski bulk. 

modified gravity. In the GBIG model we have included the Gauss-Bonnet term in order 
to modify the high energy dynamics of the universe. The GB term has been shown to 
arise naturally within a string theory context. Therefore we investigated a model where 
both the induced gravity on the brane and Gauss-Bonnet terms in the bulk are present. 
Note, in the context of string theory terms of higher order than the GB term will need to 
be incorporated at even higher energies, and earlier times. 

We have shown that the GBIG model gives some intriguingly distinct phenomenology 
from that of the DGP and GB models. With a non-zero contribution from both the GB 
and IG terms we find we have three solutions (in the zero brane tension case), as shown in 
Fig. 15. II We have a solution that starts in a finite density big bang and then self accelerates 
in the future. This is the solution of most interest. The DGP and GB terms do not remove 
the initial singularity on their own. Only when both sets of terms are present that this 
feature occurs. It has been shown that a 4D heterotic string using a one loop corrected 
superstring effective action with GB terms and dilaton and modulus fields can avoid the 
initial singularity 117011 . 

There is a second GBIG solution which also starts with the finite big bang but then 
super accelerates. This solution is therefore less relevant. The third solution which starts 
with a standard big bang singularity evolves to a Minkowski universe. The GBIG model 
has finite density, pressure and temperature but it does have a curvature singularity at the 
birth of the universe. 

The GBIG model has a severe UV-IR bootstrap which forces 7 to be very small in 
order to have high enough initial redshift. The analysis of Big Bang nucleosynthesis gives 
us the constraint that 7 must be very small. This also means that observations of Big Bang 
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Figure 5.2: Generalised GBIG models with negative brane tension in an AdS bulk. 

nucleosynthesis are incapable of discriminating between the GBIG and DGP(+) models. 
In order to discern between the two models we need to consider the earlier universe where 
the GBIG model is substantially different from the DGP(+). These differences will be 
apparent in the analysis of inflation dynamics + perturbations giving rise to formation of 
structure. The problem of structure formation in the DGP model has been investigated in 
Ref. iTTTIl . Inflation with Gauss-Bonnet terms in the bulk has been considered in Ref. [59J. 
In both cases there are substantial unresolved problems to be tackled. Combing the two 
in the GBIG model will be a difficult but important future line of research. 

In chapter |4] we looked at the GBIG model in more generality by including both 
brane tension and a non zero effective cosmological constant (p. We saw that the effect 
of including brane tension was to shift the solutions (of Fig. I5.ll) along the ii axis. The 
consequence of this is that by including positive brane tension it is possible to make the 
GBIG 1-2 start at lower densities and end with greater acceleration. It is also possible to 
make the GBIG 1-2 branches inaccessible. GBIG 3 can now accelerate at late time i.e. 
the brane tension acts as an effective cosmological constant. Including a negative brane 
tension causes GBIG 1-2 to start at higher densities and end with less acceleration. If 
enough tension is added the GBIG 1 model ends in a "quiescent" singularity. GBIG 3 
now loiters at a density given by an absolute value of the brane tension. After loitering at 
this density it will either evolve to a de Sitter solution or the "quiescent" singularity. The 
time GBIG 3 loiters for is a function of the equation of state w and 7. The exact solution 
for this time requires further investigation. This loitering solution is of interest as it does 
not require the presence of a naked bulk singularity or a de Sitter bulk, as other loitering 
brane cosmologies require. 
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Figure 5.3: Another example of generalised GBIG models with negative brane tension in 
an AdS bulk. 



The effective cosmological constant modifies the solutions in a more complex way. 
GBIG 1-2 are not altered in nature only in size. GBIG 3 is the most affected. The GBIG 3 
branch is now split in two, as seen in Fig. |5.2[ GBIG 3 can now collapse after a minimum 
density is reached, as long as sufficient brane tension is present. GBIG 4 is a new solution 
which is only present for a non-zero which is below a certain threshold. For values of 
above this threshold GBIG 4 vanishes and GBIG 1 is allowed to collapse, with sufficient 
brane tension, Fig. 15.31 GBIG 4 is only accessible with negative brane tension present. 

The are a number of constraints on 7. The most important is imposed by the fact 
that A 5 < 0. Another constraint is imposed by requiring GBIG 1 to exist. In the a = 
cases we are constrained by the density of the big bang being greater than zero. With 
the presence of brane tension we can always have this density greater than zero. Now 
the constraint on 7 comes from the fact that as we increase 7, the GBIG 1 density range 
decreases. At a certain point GBIG 1 disappears and GBIG 3 ends at a "quiescent" 
singularity of the same density as the GBIG 2 big bang (this is in the = cases; if 
7^ then GBIG 2 meets up with GBIG 4). Increasing 7 further allows evolution to 
continue through this point causing GBIG 2 and 3 to become a single solution. 

We have seen that the GBIG model has many possible evolutionary histories depend- 
ing on the values of the effective cosmological constant 0, the brane tension a and most 
importantly the contribution of the IG and GB terms in 7. Not all the solutions are 
physical. Our comprehensive analysis of the dynamics is the starting point for further 
investigations into the GBIG model, including inflation and structure formation. 
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The work on the GBIG model in chapter [3] has been cited in a number of papers 
including |[72l73]|Ml75l|7SII7l|7S|7^ work in chapter H has been cited 

inlMZTlEZilEa. 



Appendix A 
Conventions 



The conventions used in this thesis are as follows: 

• Greek letters are used for 4D indices. 

• Roman letters are used for 5D indices. 

• 5D tensor are denoted by superscript (5) with Roman indices. 

• The extra dimension is denoted by index y when in a Gaussian-normal coordinate 
system. 

• The extra dimension is denoted by index z when in a Poincare coordinate system. 

• The metric signature is: 



(-,+,+,+). 



(A.1) 



The Riemann tensor is defined such that for AdS spacetime we have: 



R 



>ab 

cd ~ 



T2 W ~ <M • 



(A.2) 
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